The 


PHILOSOPHICAL 
MAGAZINE 


FIRST PUBLISHED IN 1798 


1.1 Eighth Series No. 10 October 1956 


A Journal of 


Theoretical Experimental 


and Applied Physics 


EDITOR 
PROFESSOR N. F. MOTT, M.A., D.Sc., F.R.S. 


EDITORIAL BOARD 
SIR LAWRENCE BRAGG, 0.B.E., M.C., M.A., D.Sc., F.B.S. 
SIR GEORGE THOMSON, M.A., D.Sc., FR.S. 
PROFESSOR A. M. TYNDALL, C.BE., D.Sc., F.B.S, 


PRICE £1 Os. Od. 
Annual Subscription £10 10s. 0d. payable in advance 


ALERE 3% FLAMM@A™é. 


is, es ang-Published by 
TAYLOR & exes LTD. 


mae 


wy 
RED LION COURY-¥ 


“fy 
ey 


9 E.C.4 


The Scientific Work 
of René Descartes 


(1596—1650) 


By 
J. F. SCOTT, BA. MSc., Ph.D. 


With a foreword by H. W. TURNBULL, .a,, F.rR.s, 


This book puts the chief mathematical and physical discoveries 
of Descartes in an accessible form, and fills an outstanding 
gap upon the shelf devoted to the history of philosophy and 
science. 


There is to be found in this volume the considerable contri- 
bution that Descartes made to the physical sciences, which 
involved much accurate work in geometrical optics and its 
bearing upon the practical problem of fashioning lenses, as 
also the deeper problems of light and sight and colour. The 
careful treatment that Dr. Scott has accorded to this work of 
Descartes is welcome, is well worth reading and will be an 
asset to all libraries. Publication is recommended and 
approved by the Publication Fund Committee of the 
University of London 


212 pages, 7" x10", amply illustrated 


Price {1-0-0 net First published July 1952 


Printed & Published by 


LAY LOR |& BRAN Close 
RED LION COURT, FLEET STREET, LONDON, E.C.4 


CONTENTS OF No: I'0. 


XCI. The Self-Energy and Interaction Energy of Stacking Faults in 
Metals. By R. W. Arrren and J. S. Praskert, H. H. Wills Physical 
Maori conve University: Of Bristol, uy ayes... cece cheeceecacoh.... S85 


XCII. Proton-Excited Energy Levels in 4N. By E. J. Buras, The 
Wheatstone Laboratory, Kings College, University of London and 
D. J. Prowssz, H. H. Wills Physical Laboratory, University of Bristol. 912 


XCIII. The Heat Capacities of Chromium and Nickel. By J. A. 
Rayne and W. R. G. Kemp, Division of Physics, National Standards 
Laboratory, Commonwealth Scientific and Industrial Research Organiza- 

“Oo SURUTEIGS TCS RG Sk Ee EGGuEA Sc enn ae nr 918 


XCIV. The Drude Dispersion Formula shown to be Applicable to any 
Medium irrespective of the Polarization Field. By Sir K. 8S. KrisHnan, 
F.R.S., and 8. K. Roy, National Physical Laboratory of India, New Delhi 926 


XCV. Some Observations on the Solid State of Argon. By D. 
STANSFIELD, H. H. Wills Physical Laboratory, University of Bristol Mee ae 934 


XCVI. Density Changes in Neutron Irradiated Quartz. By P. G. 
KLEMENS, Division of Physies, Commonwealth Scientific and Industrial 
foe ares OTeTMe atom, OVONEY 00.0... 2 take yok eG e sabes deans 938 


XCVII. The Ratio of K-Capture of Positron Emission in Fluorine 18. 
By R. W. P. Drever, A. MousK and J. Scospiz, Department of Natural 


Philosophy, HO te SN We Ms Aend a SONDN Seehd ias sels > vieie ede 1e a als Gd. c/o anh Vide oe 942 
XCVIII. The Reaction Be («, n) #C. By D. B. Jamss, G. A. JonEs 
and D. H. Wiikrinson, Cavendish Laboratory, Cambridge.............. 949 


XCIX. The Cold Work Introduced during the Fatigue of Copper. 
By R. D. McCammon and H. M. RosEnsera, The Clarendon Laboratory, 
COMSTOROL Sos 6 docc0 6.6 0. 6ks:e0)-6 o.GROTELOLOND C1 DEER RC oa Sele ae 964. 


C. Correspondence : 
Periodicity in Ice-Flowers. By I. J. VAN HeeRvEN and D. J. 
Prowse, H. H. Wills Physics Laboratory, University of Bristol 967 
Long Range Nuclear Interaction in Paramagnetic Resonance. 
By J. E. BENNETT and D. J. E. Ingram, University of 
ear ait ea ie tie ese 8 Re ek 970 
The Bombardment of Osmium by Nitrogen Ions. By G. W. 
GREENLEES and L. G. Kuo, Physics Department, University 


Fes CTA Oe ee II ce hws sso ge ecto a iaiieoins 973 
Energy Levels of “A. By F. C. Frank, H. H. Wills Physical 
‘gy oetgey OR SE POLE Wee 6's cca a tem Oe aos Oates recarae 975 
IMEC WiG WarGil DOO Kemer cine 5 Soe ys, occa giene Fs o meee tp rect woes 977 


id, 
“enti cations for the Philosophical Magazine should be addressed, post-pai 
Be the Waltons c/o’ Messrs. Tayror AND Francis, Lrp., Red Lion Court, 


Fleet Street, London, England. 


‘ 
4 
ari 
" 
} 
« +4, 4 
f al ' ‘ahs 
‘ 
+ 
’ £ 
’ F hi 
, ‘ : ae th 
i i 
— 
7 
J ? Un 
‘4 i 
. r oy } a 
§ t y ‘ 
' 
r f ‘ 
. 
ah ‘ 
« * é a: 
Vix bah aot k aere oe w ¥ } 
ieee 
, 
"@ : if 
4 We 
‘ i ’ Pb id Cm 
“ 4 
i } ) while 
ris _) 
‘ 
vs 7 
¥, = , 
A 
ee ire ; Hd ‘ 
i ifs ‘ im ty P "Sat ats 
a i Pent 4 
- 4 x 
4 * 
/ ie: ce ly § C8 o iv 
; maf j ed } 
Get . ‘ . > 
} 
wit, , 4 ee ae ee “fc ba ol ie isa 
be eth ‘7 Ais 
fe ‘ Are td Ted abd “UTS danse 
; lh ibs 7 ir =) Saul 
= * ! . 7 
43 ohy ow a a 
a. i | 
7 , ~ bi k ry \ 7 5 
a j 
/ é ap: aes ‘4 
hy ae, pS RP oe Fe 
7 Y ¥ " 
7 ‘ ‘ : * ‘ vod fy tak Lae : ate 
See 
‘ by . 
; Plage a ee 
a) " $ Ped 
hi ? , np ee ye Le 
' ’ c " 
F 
. ] _ : 4 “ : 
. it AES PES) | lial Wael a 
| ¥ 
feeb hy LE ed ie ae ae 
Sin ‘ 
a 
2 
id ’ Tier Ps } i aed 
* : id | 
a ‘ * ‘ 
ad 
- 
‘A . > ' a ; 
J ; : 
. - . i 
Ls si 
. * S 4 
i 
w 
fit ¥ > 
J 5 ¥ 
; o 5 ip p ave 
“4, a ‘ ad a ke oe 
a a i A . ) 
’ on : 
t q a a 
a y 
F ’ 
13 <q 
‘ 
~~ ng 
a ; 
' . 


Pere tee 


"7 view 


ne 


[ 885 ] 


XCI. The Self-Energy and Interaction Energy of Stacking Faults 
in Metals 


By R. W. ArrreEt} and J. 8. Praskerrt 
H. H. Wills Physical Laboratory, University of Bristol§ 


[Received August 1955; and in revised form January 10, 1956] 


ABSTRACT 


A model is considered of a crystal lattice having slip of less than a 
lattice vector on two close packed planes. The effect of this slip on the 
eigenvalues of one electron wave functions extending throughout the 
crystal is calculated, using first order perturbation theory to give the 
wave functions in a perfect crystal. The self-energy of a fault in a 
monovalent metal is found to be about 20 ergs/em?. The interaction 
energy is found to be about 0-4 ergs/em? when the faults are one lattice 
plane apart and negligible when they are further apart. An unsuccessful 
attempt is made to calculate the effect of electrons with wave vectors 
near zone faces. The theory involved a clarification of the Bloch— 
Floquet theory and the use of a contour integral to sum the zeros of an 
exponential sum. 


§ 1. INTRODUCTION 


In 1950 Barrett reported that few stacking faults were to be found in 
cold worked aluminium, but that many faults could be found in certain 
copper-silicon alloys (cf. review article by Barrett in Shockley 1952, 
p. 112). Mott suggested at that time that the difference should be 
explained in terms of a higher surface energy for faults in Al than in Cu, 
this difference in surface energy being due in its turn to strong Bragg 
reflections of electrons at faults in the case of Al and only slight reflection 
in Cu (because there are no electrons near zone faces in Cu). 

Of course no measurements have been made on the energies of stacking 
faults, but there have recently been measurements of the energies of (111) 
twins in Cu and Al by Fullman (1951) (and quoted by Fisher and Dunn 
in Shockley 1952, p. 336). Fullman finds that the ratio of twin energy 
to grain boundary energy in Cu is 0-045 while in Al it is 0-21. In the 
first case the average grain boundary energy is known to be 550 ergs/cm?, 
giving 25 ergs/cm? for the surface energy of a twin in Cu. Presumably 
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the energy is about five times this in Al. Mott’s explanation of the 
difference in surface energy is expected to be as valid for twins as for 
stacking faults. ; 

This paper fills in the details of Mott’s explanation for surface energies 
of faults, with some precision for monovalent metals and very imprecisely 
for aluminium. The results are as expected. 

Of course, the experimental determination of the surface energy of a 
stacking fault is difficult so that it would be nice to know the relation 
between the energy of a fault and of a twin. On the basis of Lennard— 
Jones forces between atoms it is found that the energy of a fault is 
roughly twice the energy of a twin. On the same basis Heidenreich and 
Shockley (1948) deduce that the energy of a fault in cobalt is 7 ergs/em? 
using the known difference in energy between hexagonal and cubic 
cobalt (and forgetting that this has the wrong sign). The reasonableness 
of this figure should not be taken as an indication that the energy of a 
fault is really twice that of a twin, for the same theory (Lennard—Jones 
forces) shows that the energy of two faults separated by one layer 
(ABC|B|ABC Frank’s fault of the second kind (1951)) is roughly the same 
as a single fault. Whereas our calculation shows that two faults have 
almost exactly twice the energy of one, even though they are only one 
layer apart. 

Unfortunately, the energy of a twin is a more difficult calculation, and 
we have not attempted it. We have only calculated the energy of two 
faults M layers apart. It would clearly be possible to extend the calcula- 
tion to more than two faults, but not many more. The energy difference 
between hexagonal and cubic would again be easy. 


A Brief Summary of §§ 2 to 5 Follows 


§ 2: An electron in a Bloch function, which has a real wave vector k, 
is unfortunately not entirely scattered into other Bloch functions. Some 
of the scattered waves look like Bloch functions but have a complex k. 
The proper development of these functions is through Floquet’s theory. 
We attempt to explain this theory, its differences from Bloch’s theory, 
and the method of using it when there is one direction in the cyrstal 
that is not completely periodic. 

§ 3: It is not enough merely to consider scattering of electrons at the 
fault surface. An energy difference is required and this can only be 
obtained by working out eigenvalues that are forced on the problem by 
suitable boundary conditions. The surfaces of the crystal are therefore 
introduced and a determinant for the eigenvalues obtained. The 
individual eigenvalues do very peculiar things when a fault is introduced 
into the cystal, but by a contour integration we are able to calculate the 
average behaviour of a great many eigenvaules. 

§ 4: The general theory so far developed is applied to a metal in which 
it is supposed that the wave functions are given by the ‘ nearly free 
electron approximation ’ (generalized for complex k), i.e., first order 
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non-degenerate perturbation theory. This confines our attention to 
monovalent metals. The theory for two faults is completely worked out 
in terms of the Fourier coefficients of the potential of the perfect crystal. 
The energy change is of the second order in these coefficients and the 
interation energy decreases exponentially with M. The special potential 
Ze* exp (—q'r)/r is used round each atom, some approximations are made, 
and an easily understood answer is found. 

§5: Some early work (1950) on aluminium is briefly reviewed and the 
difficulties of obtaining any answer are discussed. 


§ 2. SOLUTIONS OF THE SCHRODINGER EQUATION IN WHICH THE 
POTENTIAL IS PERIODIC ONLY IN Limtrep Reaions 


A crystal with one or two parallel stacking faults in it, only has periodi- 
city in planes parallel to the faults. The periodicity is interrupted in a 
direction that goes through a fault. Thus the potential in the Schrodinger 
equation 
ie 
2m 


V2b-+- V(r)b= Hb(r) ee ee a ae) 


has the property 
V(r+7,)=V(r+7,)= V(r) TOrsa lr Meee ee (2.2) 
where 7, and 7, are two fundamental periods that are parallel to the 
plane of the fault. Owing to the limited periodicity in any third direction, 
it is only possible to apply Bloch’s theorem to the 7, and 7, directions. 
In order to make this possible we apply periodic boundary conditions in 
these directions 
&(r+N47,)=¢(r+N.7.)=(r) forallr. . ». °(2:3) 
A simplification of the statement of Bloch’s theorem and also of the 
subsequent mathematics is obtained if we introduce a unit vector i not 
in the 7,, 7, plane and call p the projection of r into the 7,, 7, plane 
parallel to i, so that 
rape ame eee. .  (2.4) 


Bloch’s theorem now shows that any solution of (2.1) in which the potential 
has the property (2.2), that satisfies the boundary condition (2.3) can be 
written as a sum of solutions each of which has the form 

CHDMUaRON Oe ene) ee. (2.5) 
where u(p+7;, 7)=u(p, #) and « must satisfy the equation «. N j7,=2n,7 
(cf. e.g. Shockley 1950, pp. 399-404). The periodic function u(p, x) 
depends on « but we omit explicit reference to this fact in (2.5). It is 
possible to confine « to a definite plane defined by the vectors «, and a» 
that satisfy the relations 
he RO ite ER OS OY) 


a 


Clearly «, and «, are subject to considerable ambiguity for we could add 
to them any vectors perpendicular to both Ty and 7». However, having 


302 
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fixed them in some way, we can write 

K=1N1%1/N +N o%o/N o 14; %. integers... 7. sta ae 
n, and n, are the quantum numbers of the Bloch function. It makes 
no difference how we choose the «;—the solution 2.5) is only determined 
by the integers n, and 1». 

We must now determine how w(p, 2) varies with x By hypotheses 
we can expand V(r)=V(p, 7) and u(p, x) in Fourier series with respect 
to p, 1. 

Vip, 7) =2'V (a) exp (—10 <p): u(p, v)=LA,(x) exp (—ta.p) (2.8) . 
where the coefficients are functions of x and the sum is over all « of the 
form m%1+m%_ (m, and m, integers and «, and «, can be the same 
«’s as were previously chosen—the particular choice is again of no 
importance, since only the invariant « . p occurs in (2.8). We now insert 
the expansion (2.8) in the eqn. (2.1) and assume that term by term 
differentiation is permissible. The A,(#) must satisfy a doubly infinite 
set of second order simultaneous equations ; when the « plane is chosen 
to be parallel to the p plane these equations are 

2 
—sec? 05 [Aq +2444) (e—a) . i+ A, {cos? O(«—a)?+ [(k—«) . i]? }] 
a Vi ee a . ag vis lee era 
where @ is the angle between i and the normal to the p plane. 

Now let us suppose that V(p, x) is periodic in x for some range of 2, 
e.g. between faults, then 

V —7e+r)=V,(x), for all 77 ss eee en 
Then the A,(x) satisfy a set of ordinary differential equations with 
periodic coefficients. In our subsequent work we require a complete set 
of independent solutions to these equations. The Bloch theory is of 
no use since it only gives the form of the solutions that are subject to 
periodic boundary conditions, but the Floquet theory can be applied. 
Strictly speaking the Floquet theory only applies to a finite system of 
simultaneous equations but that it should still be true for our infinite 
system (2.9) seems very reasonable physically in that it should make 
only a small error if we cut off the expansion (2.8) for |x| beyond some 
large value. 

The Floquet theory (see e.g. Ince 1927, §15.7. Notice that there is 
no need to assume that the coefficients of the differential equations are 
analytic functions of x, ie. V(r) need not be analytic ; furthermore the 
theory can obviously be generalized to a finite set of simultaneous 
equations) states that in general a complete set of solutions of (2.9) 
ean be found, all of which have the form 


A,(x)=exp (jp2)v,(x), aS Se ue 
where v,(v-++7)=v,(v). The y’s for these different solutions are the roots 
of a ‘ secular ’ equation of the form 

det [a,,(#)—exp (tu7)dj,J=0, . . , . . (2.12) 
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If the roots of (2.12) are distinct, i.e. all the exp (iu7) are different, then 
a complete set of solutions is of the form (2.11). If, on the other hand, 
there is a repeated root, then there will only, in general, be one solution 
of the form (2.11) for that root, the other solutions being of a more com- 
plicated form. This, however, does not affect our problem since we may 
take limits of the solutions as H tends to an awkward value for which 
there are coincident roots. 

There is no difficulty in showing that the a;, are integral functions of Z, 
1.e. expressible as power series in # with an infinite radius of convergence 
(see Ince 1927, §3.31 and Kramers 1935). All the y’s are therefore the 
various branches of an analytic function of #; the branches that differ 
by 27/7 clearly do not give different Floquet functions. 

We now give some of the properties of the ’s, which we shall temporarily 
write 4; These properties are essential for our later work. 


2.1. Properties of the real p; 

If one py, is real then # must be real. This can be proved by inserting 
(2.11) in (2.9), multiplying through by A,* summing over «, and inte- 
grating with respect to x from 0 to 7. 

If » ,(£’)=p’ is real and p,(L’)~p’ for any 14k then du ,(H’)/dH 0. 
Suppose that the first derivative that was not zero was the nth then 

M(H )=p' + (B-B')"A,,+. sey? 
but then 
BH=E' +A, "(upp)"... 
(Copson 1935, §6.22, p. 121). We could therefore choose a complex H 
that gave a real 1 ,(#) unless n=1. 

If »,(Z’)=p’ is real and »,(E')A~p’ for any ik then pw ,(Z) is real for 

real H in a neighbourhood of H’. Since 
,(B)=p'+(#—E')A,+(H—EH'yAot... 


we have 


B=K'+ (py By +(e rb Bet. - - « 
It follows that B,, B,,... must be real, since if we choose ,(/) real H 
must be real. The reality of A,, A,,... follows from that of B,, B,,... 
and proves the statement. 


2.2. Properties of the complex j; 


If the potential V(r) has a centre of inversion then the complex pe(H), 
for real H, are present in conjugate pairs. It is easy to show that if 
exp (iu;2)v,,(x) is a solution of (2.9) then so is [exp (—tp ,2)v,,(—x)]* so 
proving the assertion. 

These properties suggest the following classification of the He: 
(a) complex p’s with positive imaginary part and real y’s with Rosie 
derivative, (b) complex p’s with negative imaginary part and real p’s 
with negative derivative. A yu of class (a) herinafter called w+, cannot 
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possibly become a yu of class (b) hereinafter called w~, and vice versa, except 
at a branch point ; in the neighbourhood of such a point 
u(E)=p'+A,(E—E')V?+A(H—E')P/? +... 

with p aninteger. Itis easy to show that the arg A;=n7/p (n integral) 
and that therefore the p different branches of «(#) occur in conjugate 
pairs both for H<E’ and for H>EH’. If p is even there are two real p's 
(one a wt and one a p~) for H<H'(>E’) and as # increases (decreases) 
they approach p’ and turn into two complex p’s (one a »* and one a p-) 
for H>E' (<H’). If p is odd there is only one real yw (either «* or p~) 
and as £ goes through ZH’ this real » goes through yw’ without changing 
its classification. 

There are a ‘ doubly infinite’ set of u*+’s and y~’s which we can label 
by the generic symbols «, «’, ... that mean the sets 


Mh +MohKo, My'x4+My'%o,... M,;,m,,... IMtegers. (2.13) 
All out independent solutions of (2.9) can therefore be written 


(a) exp (tpg TH Vevey *(X), (d) exp (toy @ Pace (2), (2.14) 
where different values of «’ give different solutions. 

Incidentially, the group velocity of a wave that is a lmear combination 
of Floquet functions (with fixed «) with values of » that are centred 
round [iy is ih(du/dH)~1 evaluated at upp. So that ‘ut waves’ move from 
left to right while ‘.~ waves ’ move from right to left if 1 is real. If p is 
complex, the ‘y+ waves’ decrease exponentially towards the right 
while the ‘u~ waves’ decrease exponentially towards the left. This, 
then, is the physical significance of the + classification. 


§ 3. THe DETERMINANT FOR THE ENERGY LEVELS OF 
Two Stacking Faunts 


In this section we shall apply the results obtained in § 2 from Bloch’s 
and Floquet’s theorems to the physical situation of two parallel faults 
in a face-centred cubic crystal. We shall eventually obtain an equation 
for the eigenvalues of this problem in the form of an infinite determinant, 
and we shall show how this determinant may be simplified—although the 
new determinant is still infinite—and yet still lead to a determination, 
not of individual eigenvalues, but of sums of eigenvalues in a given 
range of energy. is 

First of all, however, we must specify the vectors T1, To, I, %y, and ao. 
which were left so arbitrary in § 2. It is immediately clear how we 
should choose the first three of these vectors : 7, and 7, must obviously 
lie in the plane of the stacking fault and i in the direction of the third 
period vector of the face-centred cubic lattice. In other words, 7}, Ta, 
and ir must be three fundamental periods for this lattice. If we assume 
that the faulted planes are (111) planes of the crystal, as we shall from 
now on (for reasons that will become clear from an examination of 
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eqn. (3.9) and the related discussion), then in a certain set of rectangular 
coordinates 
110 101 i 
™=— (110), 9 72= Teh irs Tq (Ol), (3.1) 
where 7/2 is the length of an edge of the fundamental cube. « , and a,» 
must now be (in order to satisfy (2.6)) 
7/2 = 


a/2 e 
(111)-+-a,(T11)],  «g=—Y[(1T1)—a,(T11)], (3.2) 


a aR 


B 


= 


where a, and a, are arbitrary. If we take a, and a, zero, the «, and nie 
are just two of the ordinary reciprocal vectors of the face-centred cubic 
lattice. If, on the other hand, we choose a, and a equal to 4, then «, 
and «, lie in the 7,7, plane, and they are, explicitly, 
Q9A/ 2 = 2H4y/ 2. = 

= COM < (AA ee EN 
These vectors make an angle of are cos (—$) with each other, that is to 
say, they are the obvious reciprocal vectors of 7, and 7, in the 7,7, plane 
and form a triangular net reciprocal to the triangular net formed by 
7, and rz. Both these choices of «,«, are useful to us. The first gives 
us an obvious correlation between xu of § 2 and the k vector of the three 
dimensional Bloch theorem for the eqn. (2.1), where we have 


V(r+ir)= V(r) 
besides (2.2): The Bloch—Floquet function that we found in § 2 defined 
by eqns. (2.5), (2.8), and (2.11) is 


Ol 


exp (tc . p+pux)Lv, (x) exp (—t« . p) Py tq 504) 
and is to be compared with exp (7k . r)u(r), so that 
k.ptpw=k.r=k.(p+ixv) . ... . (3.5) 


(from (2.4)) for all p anda. If « 1s chosen perpendicular to i(a,=a,=0), 
we can write the left hand side of (2.5) as k . r with 


k=x«-+p(111)/1/2. Merwe aes. (5.0) 
If, on the other hand, a,=a,=} then the relation between x, w and k 
is not so simple, but the mathematics of § 4, for example, is easier. So we 
shall still leave the choice of a,, a, in (3.2) until a choice is necessary. 
Since we have used Floquet’s theorem for the i direction, we can consider 
erystals in which periodicity in the i direction does not persist indefinitely. 
Let us then lay down L layers of face-centred cubic material with the 
fundamental translation periods given by (3.1). We fix the origin of 
the x coordinate so that the lattice planes are given by 


=(m+4)r or r=p+(m+4)ir, m=—L, —L+1,...—2,—1. (3.7) 
We shall only consider lattices with one atom per unit cell, and we shall 
take the nuclei to be at the positions 


r=M47,+Motet(m+4)iz, mym, integers. . . . (3.8) 
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The potential inside this block will be taken as periodic so that 

Vip+r;, 2)=V(p, 2)=V(p,#-+7), —Lr<2—— 7. > (8.9) 
This is the crucial approximation that enables us to solve the problem: 
we have assumed that the potential right up to «=0 is unaffected by 


what happens for positive x (and similarly to the left of x —L). The 
mathematical statement (3.9) carries the following physical implications : 


(a) The nuclei remain in their normal positions given by (3.8) despite 
the disturbing influence of whatever is present for x>0 (even if this is 
nothing). 

(b) The potential V(p, a) contains the potential of the valence electrons, 
so the charge density of these electrons is supposed to be unchanged by the 
presence of the boundary. This precludes the possibility of a self- 
consistent solution. 


(c) The potential due to any ions we may place in the region x>0 is 
not supposed to extend into the region «<0. 

If we were to attempt to take account of any of these points, we would 
be forced to solve a ‘ boundary layer ’ problem rather than a ‘ boundary 
surface ’ problem, a practical impossibility. 

Of cource all these difficulties are at their smallest when we choose 
to make the plane 7,, 7, the closest packed plane. This is the reason 
for choosing to consider faults on (111) planes. 

Suppose we grow on our substrate (~<0) more of the same crystal 
displaced by an amount t (lying in the 7,7, plane) with respect to the 
substrate. That is to say, we grow M layers of atoms whose nuclei 
are at the positions 


r=M47,+MoT.+t+(m+4)i7r, m=0,1,...M—1 (3.10) 
so producing what, if we may extend Frank’s (1951) definition a little, 


we may call a translation twin. If t happens to be an integer multiple 
of 7, and 7., the crystal is perfect through ~=0. If however, 


t=(7j-+7,)/3; .. © ge eo ee 
then the crystal is not perfect but is still close packed. This is most 
easily seen by projecting the position of the nuclei (3.8) and (3.10) on the 
Ty, Tg plane (in the direction (111) not i) 


[m+ (m+4$)h]7,+[mo+(m+4)4]7,, m=—L,..., —2, —1, 
(3.8’) 
[my+(m+8)3]7,+[m.+(m+3)4]7,, m=0,1,..., M. (3.10’) 


Clearly we have just removed one plane of atoms from the perfect crystal 
actually m=0 in (3.8’) or m=—1 in (3.10’). Such a erystal is said to 
have a stacking fault of the first kind in it. There are clearly only three 
types of layers, on projection on the 7,, 7, plane; In (3.8’) we can call 
them m=—3, —2, and —1 or A, B, and C. Then m=0, 1, 2,3... in 
(3:10') ara BO OeAw hae There is no translation, other than that 
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given by (3.11), which leads to a close packed arrangement different from 
the perfect crystal. The potential for «>0 clearly has the periodicity it 


V(p,%+7)=V(p,x) for 0<a<(M—1)r, eR oct) 
and is related to the potential for «<0 by 
: Wipe taaat j=) (p.a) for. —7<7=0.- 5... (3.13) 
In other words (from (2.8)) 
Ver) exp (a. t)=V (a) for —r<r<0. m(3.14) 


Furthermore, we see that, if A,(z)=exp (jux)v,(x) is a solution of the 
differential eqn. (2.9) for <0, then exp (iux)v,(x) exp (ix . t) is a solution 
of the equation for ~>0 where v,(xz) is then defined as being periodic 
from —Lz7 to (M—1)r. Of course, these two solutions of the differential 
equation, one valid for x<0 and the other for x>0, do not fit together 
at x~=0. 

On the plane z= M7 we now introduce another translation twin and 
continue the crystal to Nv. So that to (3.9) and (3.12) we must now add 


V(p, x+7)=V(p, 2) IM =e) ss. 1 (8215) 
and to (3.13) we must add 
V(p+t’, x+7)=V(p, x) (M—l)rae=—Mr, > 2. ° (3.16) 


We do not write t’=t because we wish to regard t and t’ as variables. 
This will enable us to follow the variation of the eigenvalues as slip 
occurs and also to put t’=0 and so obtain the energy of a single fault. 

In the same way as before, if exp ({uxv)v,(~) is a solution of the differ- 
ential eqn. (2.9) for «<0, then exp (iux)v,(~) exp [ix . (t-++t’)] is a solution 
of the equation for «>Mz, where v,(x) is now periodic for all « from 
—L7 to (N—1)7r if the crystal extends from —L7 to Nr> Mr. 

If we make M—=1, then we obtain a single stacking fault of the second 
kind, in which the sequence of the lattice planes is .. ABC|B|ABCA 
In other words, an additional plane of atoms has been inserted into the 
perfect lattice. From the calculation of the energy as a function of the 
integer M we shall be able to deduce the interaction energy of faults 
of the first kind and in particular the self-energy of a fault of the second 
kind (ef. the similar results for a simple solid given in § 1). 

In each of the three regions of the crystal—to the left of x=0, between 
a=0 and Mr, and to the right of z= _M@7—the solution of the Schrodinger 
equation (with given #) can be expanded in terms of the fundamental 
Floquet functions (2.5), (2.8), and (2.11) : 

Zt exp (Upbgy Toy * (#) + Ay exp (ip &)Veyqy (22) 
‘ ‘ —Lr<x<0 
2b yt OXP (ply tH Vyqrt (x) exp (ta. t)-- 
ae 4 Sby EXP (pay @)Vaqe (ee) exp (iat . t) 
| fi 0O<a< Mr, 
(3.17) 
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and similarly with c,+ and c,~ for the coefficients in the region 
Mr<a<Nvr. Inorder that the A,(«) defined in this way for all relevant x 
should be a solution of the Schrédinger equation for all «, the coefficients 
must be chosen so that A,(w) and A,(x) are continuous functions of x 
at «—0 and a=Mr. Since v,,(x) is a periodic function of x for all x, 
we have 0 ,(0)=V,(Mr), and we shall write this as just 0... Similarly 
we shall write dv,,(0)/da=v,.. Furthermore, we shall omit the 
+ and — when + and — expressions occur symmetrically. The equations 


of continuity at v=0 are 

LA Vent = 210 yi Vouyt OXP (20 - €) 
and (3.18) 
ZA (Yl Ya’ + Vo’) = 2 Oe (tpg? Vga! + Vou") exp (0 : t). 


, 


4 
Similar equations are obtained for z=Mzr. We note here for future 
reference that jy, Vyy, and Vv, are analytic functions of H and also of 
course functions of the discrete variable « that can have VN, different 
possible values given by (2.7). 

If we had one Floquet wave ‘ incident ’ on the left of x=0; i.e., only 
one a,-+ not equal to zero, then the eqns. (3.18) would suffice to determine 
uniquely the a, and the b,,*, i.e. reflected and transmitted waves (it 
being assumed that b,- =0). The eqns. (3.18) are in fact the mathematical 
expression of the so-called ‘dynamic theory’ of electron diffraction. 
These equations contain the rule for finding the wave vectors of the 
reflected or transmitted waves: If exp i[k—a-+p,t(#)(111)/1/2].r is 
the incident wave, then exp i[«—«’+p, (E£)(111)/\/2]. r is a possible 
reflected wave, where «’ can take any value of the form (2.13), but H 
is the same as for the incident wave (see (3.6) and discussion of eqn. (4.4)). 
The amplitude of this reflected wave is |a,—|? and will not in general be 
zero for any «’. This gives all the reflected waves that would be expected 
from a two-dimensional ‘cross grating’ with fundamental periods 7, 
and r,. Of course, there will only be a finite number of «”s that make 
Ly (H) real (the more the bigger £ is), viz. give undamped reflected waves. 

We are not directly interested in the waves scattered by a fault, as we 
would be if we wanted to calculate the electrical resistance due to it ; 
rather we want the allowed energy levels of the electrons. That is to say, 
we must now impose some boundary conditions at x=—L7 and Nr. 
If we imposed periodic boundary conditions, we would introduce the 
interaction energy between a periodic array of faults separated by a 
distance (L+N)r, in order to avoid introducing all these spurious inter- 
actions, we impose zero boundary conditions. This certainly corresponds 
more precisely to the actual physics of a bounded crystal and is, in any 
event, no more difficult to handle. Of course there is an interaction 
energy between a fault and the surface of the crystal; however, the 
mathematics shows that, if we make L and N tend to infinity, the sum 
of the changes in the eigenvalues due to introducing the two faults tends 
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to a definte limit. This limit is presumably the energy of two ‘ isolated ’ 
faults a distance M7r apart. 


The boundary conditions imply 
2Ay Veo EXP (—iptyLT)=0 and Zeyvyy exp [ia . (t-+t’)] 


a XO 

exp (yw ,Nr)=0. Fe er) 
The eqn. (3.18) (and the corresponding ones for c= Mr) and (3.19) are a 
set of homogeneous equations for the coefficients a,b, and c that only 
have a non-trivial solution if the determinant of the equations is zero. 
The elements of this determinant are functions of H, and so this gives us 
an equation for the eigenvalues of E (a different one for each value of x). 
The problem has two simple quantum numbers « corresponding to the 
perfect periodicity in the 7,, 7, planes, but the usual third Bloch quantum 
number (eqn. (3.6)) is no longer a constant of the motion, and our eigen- 
values cannot be classified by ». The determinant is 


— twa’ Lr 
Ue ae 0 0 
ia. (t+t’) niuo’ Nr 
0) 0) Voor © e 
ia.t 
Vee’ Veet © 0 3 90) 
ia .t (3. ) 
Wao! Wey © 0 
iam. t nia’ Mr ia. (t+t’) ,iua’ Mr 
0 Voigt © e Veo © e 


0 1D) el =€ eve Mr ia .(t+t’) elve' Mr 


Where Wyy=Yy Vo t Van" 
This appears to have six rows but only three columns, but actually 
each element should be split in two so doubling the number of columns, 


e.g. the top left hand element should really consist of the two elements 
Von? EXP (—ipytLr) and vw, exp (—tp, Lr). 

Furthermore, each element is really an cox co matrix with rows and 
columns specified by the vectors « and «’ respectively (« and «’ can 
only take discreet values (2.13)). ' 

This determinant looks somewhat less formidable if we imagine t =0 ; 
i.e. only one fault present. We can then subtract the third column 
from the second column (as written in (3.20)) and factorize the resulting 
determinant to give 


—itpwa’ Dt 
One e 0 


ta. € ptpya’ Nr 
: Cae , ia. t piva’ Mr : (3.21) 
Way’ © e 
ta .€ 
Vou’ Vow’ © 
ia.t 
Wa’ Weyer! e 


The non-zero factors exp (ia .t) and exp (i4M7r) can be removed from 
the rows and columns of the second determinant, and we see that what is 
left cannot possibly be zero for any H. It is the Wronskian of the 
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independent set of solutions of the eqns. (2.9) that we have labelled by 
an «’ anda + or —. This means that eigenvalues for only one fault 
are given by the zeros of the first determinant in (3.21); a result we 
could have obtained directly by not introducing the slip t’ on x=Mr 
in the first place. Similarly, if we put M0, subtract the last column 
from the middle column of (3.20), and factorize, we again obtain the first 
determinant of (3.21) only with t+t’ instead of t ; as we should expect. 

Returning to the general case, we see that the elements of the first two 
rows of the determinant (3.20) contain the factors exp (ju, 7) and 
exp (iu, Nr). Since we are going to make the boundaries of the crystal 
very distant from the faults, L and N will be very large (of order 10°). 
On the other hand, /—the number of lattice planes between the faults— 
will be quite small (not more than of order 10?). This largeness of L 
and N will mean that the elements in the first two rows will oscillate, or 
exponentially vary, with the energy, so producing a very dense set of 
eigenvalues. We shall, therefore, isolate these first two rows by perform- 
ing a Laplace development of the determinant in terms of them. In 
order to justify this and all the succeeding steps, we shall, have to imagine 
that the determinant (3.20) is really finite ; this we accomplish, as in § 2, 
by supposing that the v,,, are all zero for « and «’ bigger than some number. 
When we have completed all our calculations, we shall take the limit 
as this number tends to infinity. As we have said before, we have no 
mathematical justification for this step ; but we believe that it is physically 
justified to neglect the higher harmonics in the expansion of the wave 
function (2.8). It is convenient to introduce a symbol to represent the 
finite number of «’s or «’’s that we are dealing with; the symbol we 
introduce is oo. So that the determinant (3.20) has 600 rows and 
columns. To obtain the Laplace development of (3.20) in terms of the 
first 200 rows, we must select 200 columns (out of the first 200 and the 
last 20co columns of (3.20)) to make a 200 X2co determinant, which we 
then multiply by its cofactor in (3.20) and sum over all possible selections 
((400)! /(200)! of them). For example one term in the Laplace develop- 
ment of (3.20)—and a very important one, as we shall see—is 


+ a -ty TL ; 
epee) wong 0 
j — pia .(t+t’) tu N: 
f) Urn) © eta T 
=a in. t 
Vou’ Voc’ ey 0 
Oo Wey! e’@ ~t 0 
f) abe ta .t ,iuo’Mr + pia .(t+t’) atu. tM . 3.22 
‘ax’ & e Orsi ( ex z : ( ) 
0 Ww... ia. t ,ina’Mr + pia. (¢+t’ iu, +M 
ao’ © e Wry © ) eta 2 


where the middle column of the second determinant contains 2.00 columns 
while all the other columns of both determinants contain just loo 
columns. Obviously the first determinant can be factorized 


Xp it (2 pat Sng v) ; Ure 0) 
a’ a’ 


aa 
ihn (Ct) 
0 Dug 2 


(3.23) 
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Similarly the general term in the Laplace development will contain the 
factor 


exp it (2 peel-+2iei), yo ee Ree 


where the first sum means: select any co,<0ox’’s and choose in any 
way a —- or — to go with each one. Similarly with the second sum, 
with 0c0,+ 0,= 0. 

The procedure from now on depends on the number of real y,- that 
exist for the particular H we are considering. We remind the reader 
that it is possible to read off the real .,, for any given « and £ from the 
energy contours of the perfect face-centred cubic crystal. If we draw a 
line (cf. eqn. (3.6) and fig. 1) through the point k=« in reciprocal space 
that is perpendicular to the plane of the fault then each intersection of 
this line with the energy surface # gives a real ,,, for the given « and H. 
ws that differ by 27/7 are not counted as different, and the + or — is 
decided by seeing whether , increases or decreases with increasing L. 
Now it is perfectly possible for there to be ranges of H in which there are 
no real x's. These are, so to speak, energy gaps ; but notice that these 
energy gaps are dependent on x, i.e. the component of the wave vector k 
parallel to the surface of the fault. In the event that we are in an energy 
gap, so that there are no real y’s, it is immediately obvious that the 
term (3.23) in the Laplace expansion of (3.20) completely dominates all 
the other terms (3.24) for large L and NV(yu* has a positive and p~ a negative 
imaginary part). As a consequence the zeros of (3.20) are determined 
asymptotically by the zeros of (3.22). (This is a well-known theorem 
in the theory of exponential sums, cf. e.g. Langer (1931).) All this is very 
clear from the physical standpoint : any zero of (3.22) in an energy gap 
will give the eigenvalue of a ‘ bound ’ eigenfunction, i.e. a function that 
decreases exponentially on going away from the faults, but such eigen- 
values could be obtained by using only the uw, solutions in (3.17) for 
x<0 and only the »,* solutions for x> M7 and neglecting the boundary 
conditions (3.19) on z= —L7 and a=W, so giving the second determinant 
in (3.22) directly. Zeros of the first determinant in (3.22) are of no 
interest to us ; they must either be zeros of 

| Yat CXP (ty tLe) | or | yy eXP (ty Nr]. . (3.25) 
These zeros are the eigenvalues of ‘ surface states ’ on the surface x= — Lr 
or x—Nr. Mathematically it is clear that this first determinant is 
irrelevant to the calculation of the energy of the faults, because it is 
independent of t and t’. 

Now let us consider a range of H for which there are two real p's, say 
for «’=B, ie. wg are real. Then there are four terms of the form (3.24) 
that have equal modulus ; namely : 


exp i7(—pgtL+ppgtl), exp i7(—ug*L+pyp N), 
exp i7(—pg L+pptN), exp it(—vg L+pg WV), (3.26) 
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where we have neglected the common factor 

exp ir( Sy —pyth+ 2 pW) MRR AA 

a’ #B a’ #B 2 

All the remaining terms in the Laplace development can be neglected 
for the asymptotic determination of the zeros of (3.20) for large L and N. 
(The mathematical justification is as before, Langer (1931).) Similarly 
if there are four real .’s, then there are 16=2* dominant terms ; if there 
are six real p’s, there are 64—= 26 dominant terms. Thus for large L and V 
the eigenvalues are given by the zeros of a function of the following form 


A(E, t, t', N)\=2h,(H, t, t’) exp (1€,(L)N), 
l 


where the sum is over 22, 24,... values of / for 2,4,...real w’s. The 
conditions on the h, and €, required for the proof of Theorem 2 in the 
Appendix are satisfied provided we take 
E(E)=r(—pgth/N+pg-), Eq) =7(—wg-L/N +g") 
in the case of two real y.,s for example (cf. (3.26)). , is always given by 
(3.22) with the exponential part factored out. Theorem 2 shows that if 
E,(t, t’) are the zeros of 4 then 
d ? ea, 

Na ema 
where d/dt means differentiation with respect to one of the scalar 
components of t or t’ and L=O(N). In order, therefore, to evaluate 
the rate of change of the sum of the eigenvalues in (c,, ¢,) we need only 
evaluate two determinants, namely h, and h, neither of which contain 
N or L and so are rather smoothly varying functions of 2. Notice that 


Theorem 1 of the Appendix shows that the density of states is unaffected 
by the presence of the faults. 


| “log hy/hy dE 


dt &:<Ep <€. “ 


§ 4. THE ENERGY oF Two STACKING Fauuts In A MONOVALENT CRYSTAL 


In § 3 we have shown how to evaluate the change of the sum of the 
one electron energy levels when two faults are introduced into a crystal. 
The calculation depends on determining the two functions h,(#) and 
h(E), h, for example being given by (3.22). These determinants, 
h, and h, involve the wave functions of a perfect crystal (including 
those that have complex wave vectors, k) and we must therefore decide 
on what wave functions to use. In this section we shall use the wave 
functions obtained by perturbation theory, ie. regarding the periodic 
potential as a small perturbation, and we shall only go as far as the first 
order in the expansion in powers of the potential (neglect V2 compared 
with V). This is usually called the ‘ nearly free electron approximation ’ 
(cf. Mott and Jones 1936, p. 60). We shall imagine that there are 
sufficiently few electrons per atom that we do not have to. worry about 
degeneracies between different energy bands. The effect of this de- 
generacy is discussed in § 5 where we consider aluminium. 
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We have then to solve the equation (2.9) to the first order. The 
solution is of the form (2.14) where Vyq’(%) 18 periodic with period 7. 
We expand Vyv(v) and V(x) in Fourier series : 


Vea) De. exp (—niyx) 1 


Vou (#)= Yawn EXP (—niye) 
This implies that the at [V(r) of (2.1)] is 
Vio, 2) = a Vim XP [—t(a. ptnyr)) . . . (4.2) 
and the wave functions [y(r) of (2.1)] are [(2.5), (2.8), (2.11) and (4.1)] 
b(p, x)=exp [A(« « ptpyt)] 2 gory EXP [—i(%.p+nyx)]. (4.3) 
The expansion (4.1) with (2.8) and (2.9) give 


h2 
sec” 6 = {(My— Py)? —2(k—«) . i(Ugy—py) + (x—« )? cos?d 
+[(k—a) ‘ Ceara 2 V op —mY atom = EV ecg em ee (4.4) 


as a set of homogeneous simultaneous equations for v,,,,,, « and p 
specifying the equation, «’ being fixed. In this equation, H is known 
while .,, is the eigenvalue. The zero order solutions of (4.3) are plane 
waves that may be taken to be exp 2[u, «+(«—z«’) . p] where 


2nE 1/2 

Py (K—a') 2 Fb { - (Kk 2) Cog 0) s) (4.5) 

if Vop=0, ie. Vyyo=l and Vyy,—=0 if p40 and «A~a’. The wave function 
to the first order is (4.3) with vy. 9=1 and the other v,,.,.’8 given by 


h? ‘ 
Veo SEC? A me {barby —2(e—2) (Ly —py)+(«k—«)® cos? @ 


: 2H 
+[(«—a) . i]?— =} a ee (4-0) 
The part in curly brackets can easily be factorized to give 
h? 
Vow’ sec? 6 WM {iy — Mae — BY Hoe Pec —py}= Vee . (4.7) 


The expressions that occur in the determinant h,, (3.22), that we have 
to evaluate are 


Vgqi(O)=2V op 
p 
and 
Wog'(O) = tpg PaO) FY er (O) = 2V a00' p94 (Mor— DY) ee (4.9) 
p 


If «—«’ these contain a zero order term, namely | and iy, so that each 
of the sub-matrices, in (3.22) has zero order terms in the diagonal and 


first order terms elsewhere, 
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We now expand the determinant h,=A+B+C-+, where A B, and C 
are the zero, first, and second order terms, respectively. (It is simplest 
to diagonalize the matrix to the zero order first and then expand.) Then 


log hy=log A+.B/A+(C/A+B7/24%)+... . . . (4.10) 


to the second order. It turns out that A and B do not contain t or t’, 
so, for our purposes (3.28) log hy=C/A. 
After some tedious algebra we find 


log (hy/hg)=Lo +1 (t) +1, (t')+Lo(t, t'{M) = aes eee 


where J, is independent of t and t’ and 


2m\2 
I,(t)=— (=) cost O27 expili a, t) V0 Yo 


x [Kt+(a, n, m|0)+-K-(«, n, m|0)|—ee, lgeets, Loe 

2 fi 2 . ° / = Li 
LG CADE €a) cost @ & [exp (ia . t)+exp (2a . t’)—exp {ta . (t+ t’)}] 
X Ven V—o—mL-K+(a, 2, m|M)+K-(a, 2, m|M)—ce. . (4.13) 


with 
K(%, n, m|M)=exp [e(Mo—Hae M7] L(Ho” —Ho Meat Hei) 


x (Ho Ha = NY) (Ho Ha my) Pt. (40 


I,(t) gives the self energy of a fault while /,(t, t’|JZ) gives the interaction 
energy. hh, is obtained from h, by merely interchanging p,* and p,~, 
these being the only real p’s. 

The reader’s attention is drawn to two facts ; first, the energy can be 
split into self and interaction energies ; secondly, 


I,(t)+J,(t’)+-J,(t, t’]0)=J,(t+t’). 


We now imagine that all states are occupied up to some energy ZH’ so 
that we must integrate (4.12) and (4.13) with respect to E from h2x2/2m 
(Ho must be real, (4.5)) to H’; these are the <, and ec, of eqn. (3.28). 
Then we must sum over all «’s given by eqn. (2.7) inside a circle of radius 
(2mB"/h?)"?,  « and H# occur in K only, so replacing the sum over « by 
an integral, we find 


N.N. (r 
1 2 | | K dl de=L(s, n, m|M). . . (4.15) 


|u1- ao 


z [x dE= 


New variables of integration, €, yn, and ¢ are more convenient. € and 
7 are rectangular coordinates in the « plane chosen so that 


o.(h?/2mH’)1/2—(c, 0, 0), K(h?/2mH")1/2= (€, n, 0), . (4.16) 


while ¢ is chosen so that 
Hah int te)s ae ee ee) 
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The region of integration is now a hemisphere of radius 1, viz. 
Sey +o) and ¢ 0. 
The Jacobian is given by the equation 
dk dE=(4mE"?/h?)C dé dy dl. Sree aS) 
It is also convenient to introduce the following abbreviations 
a . i sec O(h?/2mH’)1/?2—f, y sec O(h?/2mH')1?2=d 


TOON OL Zi fe)\t =A; Wm ns ts 4. -~ (4,19) 
Then from (4.5) 


Mot —By =2(2mE"/h?)1/? cos 66 
by by = 20(2mE'/h?)1/? cos 0{(e— €)?— (E222) 1/2 ee (4.20) 
HoH = (2mH'/h*)™? cos O[fA6+7{(c—EP— (+) }1?], 


We can expand the square root that occurs in the last two of these 
expressions, because €?+-7?+¢?<1 and c>2 for a monovalent metal. 
The first term in the expansion is (c—€). When (c—€) occurs in the 
denominator of K (4.14) we replace it by c, but we shall not so change it 
in the exponent since this term oscillates. The integral (4.15) reduces, 
with these approximations, to 


_ N,N, 4mi” sect a fe Ht exp [(tf—c)M)] 


2mH'] c(f—nd-+1c)(f—md-+ic) 


lec alee, 4¢ 
x | | fexp LUA(ELi0)] dé dy ae. eee fe (1.51) 


The triple integral over the hemisphere of radius 1 can be worked out by 
expanding the exponential and integrating first over a semicircle perpen- 
dicular to the € direction and then over €. The answer is 


on 2 /MA\2P+1 1 
| | | exp [MMEtiL)] dé dy d= J tai & (>) p (p+2) (2p+1)" 
(4.22) 


Because we require only L++-L~ (cf. (4.12) and (4.13)) the infinite sum is 
irrelevant. It is amusing to note that this is the same answer that would 
have been obtained if (c—£) had been replaced by c in the exponential 
term of (4.14). 

The next problem is the calculation of V,,, defined by (4.1) 


Vine i { Vp, 2) dp dx / | | dp i oe (4.28) 


Again the coordinates are unnecessarily complicated, if we change to 
rectangular coordinates we introduce the Jacobian cos 0. Furthermore, 


ae Ah ee Ny, NO 
a. p+ny~= (»— aaa + ae (eptwi)y=k.r . (4.24) 
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where i’ is a unit vector perpendicular to the p plane. So our integral 
(4.23) is just 


eae i exp (ik. r) V(r) ar / far. en idee 
No nucleus is at the origin of coordinates r but there is a nucleus at 


r=7i/2 (cf. (3.8)) so we have, for a potential spherically symmetrical 
about a nucleus, 


esa | n : (—1)" 4a Ts . 
Vie ( = [vn exp (tk. r) dr= Teale rV(r) sin kr dr, (4.26) 


where Q—47r,3/3 is the volume of an atomic cell and we have used the 
Wigner-Seitz equivalent sphere. A simple potential to use is 

Vir)=Ze* exp (—g'r)ir, 7 eee 
If we replace 7, by co in the integral (4.26), this gives 
(—1)” 4nZe? 
Voen= Tikets q? ek 


from (4.24). 
Now it will be observed that J,(t) and /,(t, t’|M) (4.12) and (4.13), 


involve a summation over n and m. If we leave out all the factors 
independent of n and m the required sum is 


3 (1) {(f—nd+ie)(f—md-+ic)[o2+ (« . i—ny)? sec? +4") 
oA X [2+ (a . i—my)? sec? 0+-q’2]}-1= (h?/2mH’)? 
x CO tele) aed ed 
== (i2mE" ar |S (a; g) 2, 70 aa 


where q is defined as 


k?=o02+(%.i—ny)? sec?@ . (4.28) 


q'(W?2/2mE')2=q. . 2. 2... (4,80) 


The sum over m in S(«, 7) can be easily evaluated by a contour integral 
(see Titchmarsh (1939), §3.3, p.115). After a little manipulation we 
find 


cosh Cz cosh C'r , 
~Qas? — Qs? —i cos Fx 


sinh Cr = C sinh C'n 
Os CY Gs 
S?=sin? Fr+sinh? Or, S’=sin?F'7-+sinh? C'’n . (4.31) 
=i, R=fid, C=c/d, C'=4/(C?+Q?). 


Notice that # and C are not only dimensionless but they are also inde- 
pendent of the lattice parameter a=7,/2 and the Fermi energy LH’. 
We have, therefore, worked out d38 for seven different values of 


Q=4'7/(7/6). 


a S(a, q)=sin Fx { 
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Collecting together all our formulae: (3.46), (4.11), (4.12), (4.13), 
(4.14), (4.15), (4.21), (4.22), (4.28), and (4.29) we find that the energy 
per ne area of two stacking faults M layers apart is 

1 1 
N{Ng|t 1479] Zari 
¢ 4 2 ¢ 

as 1 or sae [d?S(a, q)}?{(2 exp [ta . t]—2) 


t=(7,+72)/3 


eee 
>| log 724 


t=0 


—exp [27(¢F —C)M](2exp [ta . t]— exp [2i« . t]—1)}+ cc, . (4,32) 


where n is the number of electron per cell, i.e. for a monovalent metal 1. 
The first part in the curly brackets gives twice the self-energy ; the 
second part gives the interaction energy. 

We have summed (4.32) over the six nearest neighbour «’s, namely 
4%, %4+%, —%1, —X%_, —%4—%. - The six next nearest neighbours 
make « . (71-+7 )/3=-+27 and so do not contribute to the energy. We 
have neglected the next nearest neighbours because d?S(«, q) decreases 
exponentially with C or ¢ or |a|. 

In table 2, we show the numerical values of the surface energy per 
unit area as calculated from (4.32) with Z=1, a=3-61 10-8 cm (the 

lattice parameter of copper) and for various values of q’7 as shown 
| (r=2-55x10-8 cm). Notice that our wave functions correspond to the 
ls band because we have used the nearly free electron approximation 
so that it would be illogical to put Z= 


Table 1. Energy of Two Faults M Layers Apart 


M 1 2 3 


Energy 2W+4W’ 4W+5W’' | 4W+6W’ 
Interaction-energy | —2W—2W’ WV! 0 


Table 2. Energy of Two Stacking Faults in Cu in ergs/cm? 


Interaction-energy 


qa Self-energy 
M=1,x102 | M=2,x105 | M=3, x 10° 
1 +20:3 +37 —28 —19-2 
2 +17-0 +33 —26 —16-7 
5 4 5-56 416-4 —15-9 — 7-45 
105; + 0-325 + 3:34 — 4-12 — 1:29 
ili a=" 0-820'<.10* + 0-498 — 0-665 — 0-178 
17-438 | + 0-94x10-6 + 0-451 — 0-604 — 0-162 
18 = 0-714 < 10% + 0-399 — 0-534 — 0-142 
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There are four points of interest in these results. 

1. The self-energy changes sign at about q’7=17-488 and for all q‘r 
bigger than this value the self-energy is negative. This is of course an 
unrealistically high q’ for it makes the potential (4.27) almost a re) function, 
at each nucleus. 

2. The interaction energy is about a hundreth of the self-energy even 
for M=1. This means that the self-energy of a stacking fault of the 
second kind ABC|B|ABC is just about twice the energy of a fault of the 
first kind ABC|BCA. For q'r=17-4377, when the self-energy is zero, 
the interaction energy is still present. 

3. The interaction energy oscillates in sign with a period of M=3 
and decreases exponentially with distance M. 

4. The self-energy for Z=1 and q'7=1 is very close to the experimental 
value of the self-energy of twins in copper (Fullman 1951). 

These results should be compared with the calculations of the energy 
on the basis of van der Waal attraction between ions, which give : 

1. A negative self-energy. 

2. An interaction energy that is equal to minus the self-energy for M=1 
and is negligible for M=2, 3,..., so that the energy of a stacking fault 
of the second kind is equal to the energy of a fault of the first kind. 

It seems wise at this point to warn the reader against an over naive 
extension of these results. It is not possible to calculate the energy of an 
arbitary number of faults by adding their self-energies and their inter- 
action energies. For example a fault on every layer gives a cubic crystal 
again and should therefore have a zero energy, not the sum of the self- 
energies. The reason for this is that the interaction energy between two 
faults given by (4.13)isnot correct when thereis another fault between them. 
It may be that the energy of a twin is nearly the same as the energy of a 
fault and that the difference in energy between N layers of hexagonal 
and WN layers of cubic is N/2 times the energy of a fault, but these results 
do not follow from this paper. 


§ 5. ALUMINIUM: THE ENERGY OF A SINGLE FAauLtT 


In this section we shall outline a method for calculating the energy 
of a single stacking fault in a polyvalent metal. For reasons that will 
become apparent, this calculation cannot be so complete as that given 
in §4 for a monovalent metal. From the very start it is necessary to put 
in numbers. We have therefore only considered aluminium. 

The wave functions, (4.3) used in §4 were based on first order pertur- 
bation theory, (4.6) (cf. Mott and Jones 1936, p. 60). For aluminium we 
must also consider electrons for which the wave vectors are close to the 
zone faces, so that, in eqn. (4.4) several v,.,,’8 are large (not just one), 
and degenerate perturbation theory must be used (Mott and Jones 1936, 
p. 61). This means that the determinant 
Ui Vy EXP (10 . t) 


h,= 
A ; 
Wey! W,4: * Oxp (ta . t) 


(5.1) 
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derived from (3.21) has large off-diagonal terms in its sub-matrices. 
We decided only to retain zero order terms, thus reducing h, to finite 
order. fy and v,,, were determined on the assumption that there were 
only two large v,., (this makes h, a 4x4). Both the real and. the 
complex solutions of the quartic equation for ly (cf. Mott and Jones 
1936, p. 60, eqn. (40)) were obtained. The solutions for LHy—only the 
real ones of course—are shown in fig. 1. We have plotted a cross section 
of the energy contours. The cross section. is taken through the centre 
of the zone and containing the vectors (I11) and (111). Note that we 
are now taking the « plane perpendicular to (011) so that the relation 
between x, 4. and k is given by (3.6). The numerical values for the 
aaa gaps at the hexagonal and square faces were taken from Leigh 


Cross section of the energy surfaces in aluminium. The full curves are for the 
upper band and the pecked curves are for the lower band. Contours are 
at intervals of 0-05 Rydbergs. A corner of the zone boundary between 
a square (020) and a hexagonal (111) face is shown. The line (111) is 
the intersection between the plane of the figure and the « plane, which 
has an (011) normal. The intersection of the line (111) starting at a 
given x with the energy surfaces gives the value of » (#). Note that for 
E>0-85 Rydbergs there are, for most «, four real y’s. While for 
E<0-85 Rydbergs there are only two real y’s. The curves have been 
slightly adjusted near the edge between the square and the hexagonal 
faces. 


A reasonably complete calculation would involve the computation of 
h, for many values of « and H over the required range (up to about 
0-9 Rydbergs). Actually we selected 3x 6x’s each one lying in the 
direction of one of the 6 zone faces. Five values of the energy, 0-70, 
0-75, ... 0-90 Rydbergs, all lying in the upper band were used for each 
of the three essentially different x’s. 
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In fig. 2 we show the results of a typcial calculation of the real and 
imaginary parts of h, for |«]=0-70 atomic units and « in the «,, ie. (111) 
direction. The six graphs are for six values of «,.t:0,1,2,...5. 
The curves for, .t=a«, . (71+7,)/8=47/3=4-2 (3.11), and «,.t=0 are 
the only ones of immediate interest, but the intermediate ones are 


Fig. 2 


ivan 
mala 


&.t = 3 


The real and imaginary parts of h,(2, «, t) (ordinates) plotted against 
E (abscissae) for a given « and for w, .t=0,1, 2,...5. 


necessary to get the correct branch of arg h,=are tan (S3h,/Rh,). It 
will be noticed that for this particular « h,=0 at the bottom of the band 
when «,.t=0 and 3-8. This zero of h, actually moves out into the 
energy gap between the two bands (for this «) and returns after 


t= (3-8/4-2)(7,+7)/3. 


This zero in the energy gap is, as we explained in § 3, the eigenvalue of a 
bound surface state that has for its wave vector parallel to the fault the 
value «—[(111) . «/3](111). We found no surface states for the values of 
« we chose when t=(7,+-7,)/3. 
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We estimated that the energy change due to the fault is about 
2600 ergs/cm? for «’s near the hexagonal faces and about —1600 ergs/cm? 
for «’s near the square faces. So giving a total energy of the fault of 
1000 ergs/em? for electrons in the upper hand. We do not have much 
faith in this calculation, because it assumes that the energy change is 
the same for all x’s in the sectors terminating at the faces. It seems very 
likely that for x’s near the edge formed by two faces the energy would be 
entirely different, i.e. somewhat between +2600 and —1600. 

Owing to the use in the total wave function, (2.8), (2.11), and ‘ey Ni, 
of complex wave vectors k in the direction perpendicular to the fault 
there will in general be a charge accumulation (positive or negative) in 
the region of the fault. Not only will this violate our condition (3.9) ; 
it will also give an electrostatic contribution to the surface energy. 
The electrostatic energy will be infinite unless the surface charge is 
screened, and the exact magnitude of the energy will depend on the 
rapidity of the screening. Suppose, for example, the charge density 
after screening was 


p(2)=polexp (—Ay2[)—(Aa/Ax) exp (—Aglz|)]_- - (5.2) 
where z is the perpendicular distance from the fault. Then the self-energy 
of this charge distribution is 

Qi=re)? 

* AAs(Ay-FAg) °° 

per unit area. If we take py»=10-?, A,;=1, and A,=4 in atomic units, 
the energy is approximately 200 ergs/em?. With py=10-?, A,;=1, and 
A,=+ the energy is greater by a factor of 6. These energies are of the 
game order of magnitude as those due to scattering of electrons by the 


fault. 


(5.3) 
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APPENDIX 
The Zeros of Exponential Sums 


In this Appendix we shall study the behaviour of the zeros of a finite 


exponential sum, say 
N 


A(H,t, N)= &X h,(H, t) exp (16, (#)N) 
t=1 
as a function of the parameter ¢ for very large positive values of N. 
The steps to the proof of Theorem 2, which is used at the end of §3, have 
been labelled Lemmas 1, 2, and 3. Lemmas | and 2 are independent of 
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one another, and give rather different bits of information on the zeros 
of 4. Lemma 3 uses both these bits of information and the two theorems 
follow directly from it. 

The following conditions are common to the whole Appendix (notice 
that no mention need be made of the parameter ¢ until Theorem 2). 


1. g (EB), h(E) and €,(2) are regular functions of # in a region D of 
the H plane which has the closed interval (¢,, ¢,) of the real axis in its 
interior. 

2. The €,(Z) are real for real values of # in D and can be labelled so 


that 
&(e) << & (6) <6 le); 1Al or 2 


for all « in the closed interval (e,, €,). We shall always write H=e«+27 
so that « and 7 are the real and imaginary parts of LZ. 


Lemma 1. With the conditions 1 and 2 


ce yay _J OxP Lies B)N Wy (E)+O(exp —anN)] 150 
ees AUT Gat ah =| exp [iég(B)NIIh,(Z)-+Olexp bnN)]} <0 


for |n|<yo>0 and «y—yy)<e<e€g+7- @and bare positive (non-zero) and 
independent of H and NV. 
We only consider the case of 7>0 


Xh, exp (t€,N)=exp (7€N)[hy+2"hy exp {i(E—§)N 3] 
but 
LE, (L)—&,(E)=LE(e)—€s()] alr (e)—€1' (2) + O(y?) 


and from condition 2 there exists an a>0 for which €,’—€,'>2a for « in 
(ey—7', €2+7’) for 7’ (>0) and small enough for all 7. It follows that we 
can find an 7’ (>0) so small that the term O(y?)<ay for 0 <<yj<y7"’ and so 


lexp {tL€,(Z)—&,(Z)|N }|<exp (—ayN), O0<97<%<0 


and €4—ny<e<ey+ where 7 is the smaller of 7’ and 7’. If yn is at 
the same time chosen so that this rectangle lies inside D then the |2,| 
are bounded and the lemma is proved. 

Corollary. The zeros of A(H, N) for large N either tend to the zeros 
of hy (17>) or ha(n<0), if |y|<m or are within a distance of O(1/N) from 
the real axis provided h, and h, are not zero on the real axis. At a point 
on the real axis where h, or h, has a zero the zeros of J lie within a 
distance of O(log N/N) of the real axis. 


Lemma 2. With the condition 1 and provided not all h, are zero at 
the point He the number of zeros of 4 within the square of side 2A/N 
(A arbitrary but independent of N) centred at « is bounded (as N + oo). 
Furthermore, if circles of radius C/N (C<0 as small as one likes but 
independent of NV) are cut out around each zero of A then in the remaining 
region of the square |4| has a positive lower bound (N big enough) that 
is independent of J. 
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Let us change the variable from E to z, H =e+2/N. Then 
A=Th,(c) exp [ié,(e)N] exp [é€/(e)z]+O(1/N),  z<A/2. 
Consider the first term and write it Ya, exp (ib,z)=f(z) with 
a=h,exp (i€,N) and b,=€,'(e). 


f(z) has a finite number of zeros inside the square of side 24 centred at 
z=0 for any given a,, i.e. given N. If we remove from the square circles 
of radius C' centred at each zero then | f(z)| has a positive lower bound in 
this region. The number of zeros of f(z) and this lower bound depend 
on NV through the a, but for all N the a, only vary over closed regions 
since |a,|=|/,|, so the number of zeros has an upper bound independent 
of N and |f(z)| has a positive lower bound independent of N. Now 
consider only those N’s that are so large that the lower bound of |[f(z)| 
dominates the term O(1/N) in 4A. Then from Rouché’s theorem 
(Titchmarsh 1939, §3.42) 4 has the same number of zeros as f(z) and 
furthermore |4| has a positive lower bound independent N in the cut out 
region. 

Lemma 3. 

29 (EZ) exp [té (H)N] a { ¢ Me Ye 
PEED RAW NI r,t | ng OO) 

where the closed contour I" is composed of the two contours I’, and I. 
I, starts at «,, ends at ¢, and is always above the real axis except at the 
ends points. J starts at ¢, ends at <,, and is always below the real axis 
except at itsend points. /,is not zero on J’, nor is h, zero on J’_ including 
the end points. h, exp (i€,N) is not zero on I. g,, h,, and €, satisfy 
conditions 1 and 2. 

Consider 

2g, exp (2€,N) —? |= poe ( J(Gihy—9 11) exp [1(E.—§1 NI ag 

ie Eee exp (#E,N) hy ST, ahh, exp [4(€,—€)N] 
By Lemma 1 the numerator is of order exp (—a7N) and, except in a 
neighbourhood of radius A/N of «, and «, (A sufficiently large), the 
denominator is dominated by h,? and so has a lower bound. In this 
neighbourhood of «, and «, ihe denominator also has a lower bound 
independent of N by Lemma 2. An upper bound to the integral on the 
right is therefore 


B | exp (—anN)|dB| 
Te 


where B is independent of N. This integral is of O(1/N). The proof for 
I_ is similar. 

Theorem 1. The number of zeros of 4 that lie inside the rectangle with 
corners €,-++in,) and cae (no>0 but sufficiently small) is 


=[- [é9'(e)—£1'(€)] de +O(1) 
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provided h, and &, satisfy conditions 1 and 2 and h, and hf, are not zero 
at €, and €y. 
The number of zeros inside the rectangle is 
1 ¢ A’/AdH 
2rt J p 

where I" is the boundary of the rectangle except where there is a zero 
on the boundary of the rectangle except where there is a zero on the 
boundary when I" goes inside the rectangle to avoid it. Since 


A’=S(h, +hpé,N) exp (1€ NV) 
and the conditions of Lemma 3 are satisfied with either g=h, or hie! 
(we also modify I" so as to avoid zeros of h, and h,) we have for the number 
of zeros 


oh hs hisy' ap nal hoé, 
2riJr, hy 2ntJ) pr 
which proves the theorem once it is noted that the integrands of the 
I, and I_ integrals have no singularities and so can be deformed to 
lie on the real axis. 

We now consider the variation of the parameter ¢ contained in the 
h, and impose the following condition in addition to 1 and 2. 

3. h,(H,t) and oh,/ot are regular functions of # in D and continuous 
functions of ¢ in some interval. 

Theorem 2. If the zeros of A inside the rectangle with corners €,-+7n9 
and. €,+inp (jo >9 but sufficiently small, cf. Lemma 1) are at the points 
L(t), p=1, 2,... for a given value of NV, then 


d Les eee 
Te Bal= rae log 7ide-+O(1/N). 


The zeros are to be counted with their correct multiplicity. Besides the 
conditions 1, 2, and 3 on h, and &,, h,(H, t) must have no zeros with 7>0 
and h, must have no.zeros for 7<0. (The result could easily be modified 
with the weaker condition that , and h, should not be zero at €, and €,.) 


i / 
2B y= xm] 24 /A dE, 
where J’ is as in Theorem 1. So 
d 1 
log 4at=5—([ Bo oO 
3 


a 


0? 
Teo a | aa 
integrating by parts we find 


(ae Aa Levine got 
We now use Lemma 3 with g,=0h,/dt and we find 
d ia d 
5 oEp=— aye log h ab oo log hy dE +0(1/N). 


The theorem follows AGS the appropriate modification of the contours. 
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ABSTRACT 
The nuclear energy levels of !4N have been investigated up to an energy 
of 7-7 mev by inelastic proton scattering. All the known levels appear to 
participate in the scattering. The energies of the levels have been 
determined as : 7:60--0-02, 7-40--0-02, 7:03-+0-02, 6-46 -+-0-02, 6-23-+0-02, 
5°83-+0-03, and 5-69-+0-03 Mev. Some evidence was obtained for a 
level at 5-95 Mev, and a probable new level was observed at 6-60=-0-04 Mev. 


Tue scattering of 9-5 Mev protons by nitrogen, accelerated by the 60 in. 
Birmingham cyclotron, has been studied by Freemantle ef al. (1954) 
using the scattering camera described by Burrows et al. (1951). Nuclear 
energy levels were detected at 2-3, 3-9, 4:9 and 5-1 Mev, as well as at higher 
values. These higher levels have now been examined more carefully 
with plates from the same exposure and also from an exposure made in 
October 1955. 

A study has been made of the factors limiting the application of nuclear 
emulsions to the detection and measurement of the energies of particles 
inelastically scattered from a target. A micrometer eyepiece (x15 
Kellner) was used which allowed readings to be made to 0-01 mm 
(=1 division) in the focal plane. In conjunction with a x 95 objective, 
it was found that repeated settings on the beginning and end of the same 
track gave the range with a mean deviation of 3 divisions, i.e. about 0-4 p. 
An error was found to arise from the determination of the position of the 
first grain. This was aggravated by a concentration of background grains 
near the surface. The origin of this is obscure and it has not been found 
possible to eliminate it. To facilitate scanning in earlier work, where 
range was of importance only in order to assign a track to a well defined 
group, the surface had been rubbed with cotton wool dipped in alcohol to 
remove this ‘deposit’. The rubbing was not so severe for the angles 
larger than 90° used in the present work, and it was allowed for by adding — 
a constant amount to the projected range of every track such that the 


ee 
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Q-values of well-defined groups agreed with well-known levels found in 
other reactions. The maximum addition was 16 divisions which is 
equivalent to the removal of about 6 divs in depth (=0-7 ). In the more 
recent exposures no rubbing took place; this correction was therefore 
not necessary, but the surface conditions made scanning very difficult. 
A study of the average distances between successive grains for 81 tracks 
on a ‘rubbed’ plate, starting at the surface at an angle of dip of 21°, 
led to the following results : Ist to 2nd grain, 10 divs ; 2nd to 3rd, 10 divs ; 
3rd to 4th, 8 divs; 4th to 5th, 6 divs. Thus an error in the determination 
of the first grain can lead to a significant broadening of the histogram 
for a given level. 

Straight attacks found within an angle of --ca. 4° of the mean direction 
were accepted, the micrometer eyepiece being aligned to each track. This 
criterion is an individual one and means that the number of tracks in the 
various groups are only a very approximate indication of the cross section 
for the corresponding energy levels. The mean direction is not the same 
as the line of traverse for a given angle, as can be easily verified from a 
consideration of the geometry. The area to be scanned for a given histo- 
gram needs to be small in both directions otherwise complicated correc- 
tions become necessary. Two effects limit the useful length of traverse, 
(1) the variation in the angle of dip, affecting the projected length, and 
(2) the variation of the correction for the stopping power of the gas in the 
camera. These effects act in opposite directions and for laboratory angles 
near 120° the latter proved to be about twice the former at a distance of 
about 5-6 cm from the point of scattermg. The optimum width of the 
traverse is determined by the variation of the energy of protons with 
scattering angle. For a given Q-value, this variation is given by 


E,=E, [a cos 0+ (a? cos? 6-+b+c¢Q/H,)!? ? 
where EH, is the energy of the incident particle of mass /,, 
4 a=M,M,/(M,+M,), 


where M, is the mass of the scattered particle (=, for simple scatter- 
ing) and M, is the mass of the recoil nucleus, 


b=(M,—M,)/(M a+ Ms) ’ M;/( M,+Ms), 


and @ is the scattering angle in the laboratory system. In the present 
work it was found that the expected variation in track length for a change 
in 0 of -L0-5° (ca. 0-5 mm traverse width) was equal to that expected for 
a traverse length of 5 mm, the sum of the two effects amounting to 
about +0-25 , or -—-2 divs. For some histograms, where large areas 
were needed, the variation of @ was about +1°. 

The width of the scattering gap introduces another factor limiting the 
resolution. A given point in the traverse is able to accept protons 
scattered from a finite length of the beam, owing to the width of the gap. 
For a distance of about 5-6 cm, a mean scattering angle of 120° and a gap 
width of 4 mm, the range of scattering angles accepted is 121-7°-118-2°. 
This will therefore introduce a further spread in the measured lengths of 
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upto +3divs. The width of the traverse is too small to make a significant 
increase in this correction. On the other hand, there is a variation in the 
correction due to the stopping power of the gas, arising from the different 
points in the beam at which scattering can occur. For the above 
conditions, this amounts to a path difference of ca.+1 mm, and a corre- 
sponding correction due to stopping power and angle of dip of -L0-6 divs. 

An important limiting condition is the straggling in the range of the 
particle. This is about +2-5% for a 400 divs track and +10% for a 
100 divs track (Rotblat 1950). So we may take it as +10 divs. The 
maximum total spread resulting from the area of the scan is +2 divs, 
and from the gap width +3-6 divs. These give a total standard deviation 
of about +11 divs. Although the maximum error due to the selection of 
the first grain can amount to +10 divs, the average value is considered 
unlikely to exceed --3 divs. Thus it can be seen that the major contribu- 
tion to the standard deviation is due to the straggling, which is inherent 
in nuclear emulsion techniques, and not due to the design of the camera or 
the errors in scanning. 

RESULTS 

Range histograms of tracks measured at four different angles, 140°, 
130°, 120° and 60° are shown in the figure. Measurements have also been 
made at 80° and 110° and they confirm the level scheme given by the 
other angles. The mean proton beam energy determined from the 
range of the elastically scattered particles was 9-45+0-01 Mev. The 
range energy relation of Gibson ef al. (1954) was used. The correction 
made for the loss of energy in the gas was based on the values of Brolley 
and Ribe (1955). 

It can be seen that levels are present at 7-60 (at 60°), 7-40 (at 60°), 
7:03, 6-46, 6-23, 5-83 and 5-69 Mev. There is also a possible level present 
at about 6-60 Mev which had been suspected previously (Prowse 1954). 
The uncertainties amount to --0-02 Mev except those on the 6-60 Mev 
level and the 5:83-5:69 Mev doublet where they are 0:04 and 0-03 mev 
respectively. These errors are derived from the statistical uncertainties 
plus the error in the beam energy. 

The calculated position of the 6-60 mev level does not coincide with the 
experimental position at all angles ; indeed, it appears at a higher range at 
low angles and at a lower range at high angles. This would apparently 
indicate that the level belongs to a lighter nucleus than nitrogen, e.g. 
carbon (the exact excitation energy would then be slightly different). 
If this were the case it would have been due to a gaseous impurity in the 
gas used in the scattering camera. Ordinary commercial nitrogen was used 
of 99% purity. Possible contaminations are carbon dioxide and oxygen. 
An oxygen impurity would result in a group of elastically scattered particles 
which, although not resolvable from the '4N group, would have the effect 
of increasing the apparent beam energy at high angles. No such effect 
was observed. Similarly a !2C impurity would lower the apparent beam 
energy at high angles. No tracks were observed which would be con- 
sistent with the 4-4 mey level of 0, The loss of energy in the gas was such 
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Range histograms of the tracks obtained at various angles. The excitation 
energies are the average values obtained from measurements at six 
different angles, except 7:60 (60° only) and 7:40 (60° and 80° only), 
and 5:95. The results have been smoothed by adding to the number of 
of tracks in a given 5 divs interval the average of the numbers in adjacent 
intervals. The scales on the ordinates are marked in multiples of 
10 tracks, 
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that no tracks would have been expected from the ey mev level of ?C which 
is known to participate in proton scattering (Hossain 1955). There remains, 
therefore, an anomaly which cannot be easily explained although it appears 
likely that an additional level at 6-60 Mev in '4N would provide a solution. 
No such level has been previously reported from any other reaction. 

A level at 5-98 Mev has been reported by Sperduto and Fader (1953) 
from the !O(d, «)!4N reaction and also possibly by Benenson (1953) 
who obtained 6-1--0-1 Mev from the @C(d,n)!4N reaction. It has not 
however given rise to a y-ray in this or in the *C(p, y)'*N reaction. In 
the figure the positions where it would be expected to appear have 
been marked and it can be seen that if it exists the cross section for proton 
scattering is low. The best values obtained for other levels compare well 
with those derived from y-ray transition energies, quoted by Azjenberg 
and Lauritsen (1955), of 5:68-0:01 and 5-83+0-02 Mev. The 6-23 Mev 
level has been observed only in the C(d,n)14N reaction (Benenson 1953) 
and the error quoted was 0-05 Mev. The present results confirm this 
value and reduce the uncertainty. This level has also given rise to a 
6-23-.0-03 Mev y-ray. ‘The same is true of the 6-44 Mev level which has 
been observed to lead to a 6-45 Mev y-ray from 1C(d,n)!4N. The 
7-03 mev level has given a y-ray of 7-05--0-04 mev although its excitation 
from neutron group measurements is 7-:00--0-04 Mey. Similarly the 
7-40 mev level has given a neutron group which results in an excitation 
energy of 7:50-+0-04Mmev but the y-ray energy was 7:30-+-0-05 mev. 
The 7:60 Mev level has previously been reported at 7-72-+0-04 mevy, 
again from the (d,n) reaction, but y-rays have not been observed. The 
levels of lower excitation, 5-10, 4:91, 3-95 and 2:31 Mev have all been 
observed in this experiment, the 4-91—5-10 Mev doublet being resolved at 
most angles despite the low cross section of the 4:91 Mev component. 
A figure of 3-0-+1-0 can be given for the ratio of the cross sections for 
scattering leading to these levels at this energy. No tracks were observed 
which could have resulted from a level at 3-15 Mev which has recently 
been reported by Quinton and Doyle (1956). In conclusion it may be 
stated that in resolving the levels at 5-69, 5-83, 6-23 and 6-46 Mev the 
limit of resolution has been reached for the scattering camera with the 
present aperture using photographic plates. (Similar resolution has been 
obtained by Burcham et al. (1953) using the same camera to investigate the 
doublets of 1°O.) Further investigations are proceeding to see whether, 
by slightly tilting the photographic plate in the microscope and thus 
rendering the proton tracks horizontal, better results might be obtained 
in judging which is the first grain of the track. 
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ABSTRACT 


Heat capacity measurements have been made below 4-2°K on extremely 
pure specimens of chromium and nickel. In both cases the measurements 
give y values lower than those obtained by former observers. The 
values of Debye temperature also differ significantly from those of previous 
work, it being found that 9=630°K for chromium and 6=440°K for nickel 
in contrast to the older figures of 418°K and 413°K respectively. These 
new results help to resolve the apparently anomalous nature of the @ 
versus 7’ curves for the two metals. The influence of the magnetic 
heat capacity on the results for nickel is also discussed. 


§ 1. INTRODUCTION 


MEASUREMENTS of the heat capacity of chromium have been made in the 
liquid helium region by Estermann et al. (1952) using a specimen of 
approximately 99-9°, purity. Their data, when analysed in the usual 
way into a linear and cubic term, give a value for the Debye temperature 
of 418°x.§ As the results of Anderson (1937) on the heat capacity of 
chromium from 56 to 300°K show that the Debye temperature increases 
from approximately 450°K|| at room temperature to 500°K|| at 56°K the 
figure obtained by Estermann et al. implies that the curve of 6 versus 7’ 
has a maximum other than at 7’=0, contrary to theoretical expectations 
(Bhatia and Horton 1955). The present work on chromium was under- 
taken to investigate this anomalous behaviour and also, by using a speci- 
men of high purity, to provide more accurate data concerning its electronic 
heat capacity, previous work (Rayne 1954) having shown that the y value 
for transition metals is quite sensitive to small amounts of impurity. 


' af es at Westinghouse Research Laboratories, East Pittsburgh, Pennsylvania, 
ft Communicated by the Authors. 

§ After submitting this paper, the authors learnt of a more recent measure- 
ment of the low temperature heat capacity of chromium (Wolcott, N. M., 1955, 
Conference de Physique des Basses Temperatures, p. 286). His value of €@=585°K 
is in good agreement with that found in the present work. 

|| These values of Debye temperature have been corrected for the electronic 
heat capacity of chromium and also for the difference between C, and Cy. 
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Keesom and Clark (1935) have measured the heat capacity of nickel, 
of 99-8% purity, from 1 to 19°x. For temperatures below 9°xK, they found 
that their data could be represented fairly well by 6=413°K, this figure 
-being derived from the room temperature elastic constants. Above 9°xK 
their results, taken in conjunction with heat capacity data extending to 
higher temperatures (e.g. Busey and Giauque 1952), indicate that the 
Debye temperature rises to a maximum of 450°K at about 18°K and then 
falls to a constant value of 390°K at high temperatures (Clusius and 
Schachinger 1952). According to Bhatia a face-centred cubic metal can 
exhibit such a maximum at 7’'=0 only if its elastic constants satisfy the 
condition 

as Cu—Cp—Cu 30) 
Cus 

The low temperature elastic constants of nickel have not been measured, 
so that it is not known whether the above condition is satisfied or not. 
For copper, which is the element next to nickel in the periodic table and 
which also has a face-centred cubic structure, recent measurements 
(Gaffney and Overton 1954) givec=—1-365in agreement with the fact that 
the @ versus 7’ curve has a maximum at 7=0. One might thus expect 
that nickel would behave in a similar fashion. The present work on 
nickel was undertaken to investigate this possibility and also to provide 
accurate data on the electronic heat capacity of pure nickel. 


§ 2. EXPERIMENTAL DETAILS 
(a) Apparatus 

The calorimeter employed in these experiments uses a mechanical heat 
switch, the possibility of extraneous effects arising from the use of 
exchange gas being thus avoided. Details of the apparatus are given 
in fig. 1. 

Temperatures in the range 1-5—4-2°K are obtained by controlled pumping 
on liquid helium in the inner chamber A, which is filled from the main 
bath by means of the needle valve B. The working space below the 
inner chamber is normally evacuated but provision is made for admitting 
exchange gas to calibrate the resistance thermometer at the conclusion of 
a run. 

Within the working space the specimen C is suspended by nylon threads 
from the disc D, which is thermally anchored to the bath by a copper 
strap E. Movement of the disc and specimen is effected by a stainless 
steel tube connected to a screw and bellows mechanism F at the top of the 
cryostat. By pressing the specimen against the flat bottom of the can, 
the sample may be brought to the inner bath temperature. Initial 
cooling from liquid nitrogen to liquid helium temperatures takes 
approximately 10 minutes. 

The heater for the specimen is a 5000 ohm manganin resistance, which 
is wound on a copper former G screwed into the top of the specimen, 
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A 10 ohm Allen-Bradley carbon resistance, mounted in a copper sleeve H 
screwing into the side of the specimen, is used as a thermometer. Electrical 
connections to both the heater and thermometer are made with fine 
manganin wires, which are wound several times around the copper rod 
J to effect thermal contact with the inner bath and hence to reduce the 
heat leak to the specimen. 


Fig. 1 
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When the disc is raised so that the specimen is hanging freely, the only 
heat leak to the specimen arises from conduction down the thread and 
electrical connections, stray radiation and vibration. By keeping the 
temperature difference between the specimen and inner bath small, the 
stray input can be maintained below 100 ergs per minute. This enables 


the average random error in the individual heat capacity determinations 
to be reduced to about 1%. 


The Heat Capacities of Chromium and Nickel 921 


The electrical and timing circuits are quite conventional and will not 

be described here. 
(b) Specimens 

A specimen of ductile chromium (Wain et al. 1954) of 99-998°, purity 
was used in this work. The specimen was prepared by the Common- 
wealth Department of Supply in the form of a cylinder weighing 
185-4 gramme. The nickel sample was a cylinder of mass 39-85 gramme, 
prepared by Johnson Matthey from spectrographically standardized 
sponge of 99-999°%, purity. 


Fig. 2 
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§ 3. RESULTS 


The experimental results, corrected for the heat capacity of the thermo- 
meter and heater and also for deviations from the 1948 temperature scale 
(Erickson and Roberts 1954), are presented in figs. 2 and 3 as plots of 
C/T versus T?. For a normal metal this plot should give a straight line 
whose ordinate at 7’=0 is the value of y for the metal and whose slope is 
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related to the Debye temperature by the relation 


12 7*R 
oare 68 : ce 


R being the molar gas constant. 

The straight lines shown are least squares fits to the points, the 
maximum deviation from the line being about 1:5%. The table gives 
the corresponding values of y and @ together with the results of previous 
observers. 


eT 


Chromium Nickel 
Observers ; | 8 7 @ 
(Joule mole) (°K) (Joule Mole) (°K) 
(deg.— « 10%) (deg.—* x 10%) 
Estermann et al. 1-59+0-037 | 418 +20 7 
Keesom and Clark 7340-2 7 | 413+15 F 
Present authors 1-40+0-01 ¢ | 630 +30 f | 7-02 +0-06 f | 441415 f 


+ Errors (random only) estimated by present authors from plots of experi- 
mental data. : 


¢ Errors (random and systematic) estimated at a 90% confidence level. 


§ 4. Discussion 


The present work on chromium shows that its electronic heat capacity 
is very low for a transition metal. This result is in qualitative agreement 
with the data of Estermann et al., who discuss the relation of the low y 
value to possible theoretical forms for the band structure of chromium. 
Even if their work is corrected for errors in the 1948 temperature scale 
(Erickson and Roberts 1954), however, there is stilla significant discrepancy 
between the two values of y ; it is possible that this could be caused by 
impurity effects. 

The new @ value for chromium is much higher than that previously 
obtained and gives a @ versus 7’ curve of the form shown in fig. 4. As the 
data of Anderson only extends to 56°K the dotted portion of this graph is 
hypothetical, but it is clear that the new result is not inconsistent with 
a curve having a maximum at 7’'=0. It is considered that the present 
# value is the more reliable figure and that the previously implied anomaly 
in the 4 versus 7’ does not exist. No adequate explanation can be 
advanced for the discrepancy between the values of Debye temperature, 
as neither impurity effects nor errors in the 1948 temperature scale could 
materially alter the older value. Since the most recent measurements on 
copper (Corak ef al. 1955) and on magnesium (Smith 1955) also give 
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similar disagreement with the results of Estermann ef al., it is possible 
that there were undetected errors in the latter’s experiment. 

Reference to the table shows that the agreement between the present y 
value for nickel and that obtained by Keesom and Clark is within the 
combined limit of error. It is to be noted, however, that the present 
determination is considerably more accurate, particularly since the 
estimated uncertainty in the older work does not include systematic 
errors. 

Fig. 4 
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Owing to its large electronic heat capacity, it is not possible to determine 
the value of Debye temperature for nickel with any great accuracy from 
measurements in the liquid helium region alone. ‘To draw any conclusion 
from the present work regarding the 0 versus 7' curve for nickel, it is 
necessary to make use of existing heat capacity data above 4-2°xK. 
Unfortunately, most of these measurements have been made on specimens 
of high or unknown impurity content and hence cannot be related to the. 
present results owing to the attendant uncertainty in the electronic heat 
capacity, which constitutes a large fraction of the total heat capacity at 
low temperatures. The only data on the heat capacity of high purity 
’ nickel above 4:2°K is that of Busey and Giauque but their work only 
extends to 13°K. The present values of y and 6 fit their low temperature 
results quite well, however, as may be seen from the fact that these para- 
meters give a calculated value of heat capacity at 15°K of 0-183 joule 
mole! deg-? as compared with Busey and Giauque’s smoothed value of 
0-180 joule mole-! deg~®. It is thus reasonable to conclude that the 
value of 6 has a constant value of 440-+-5°k up to 15°K. The resulting 
@versus 7’ curve then has the form shown in fig. 4, the anomalous maximum 
found by Keesom and Clark being no longer present. After correcting the 
latter’s data for errors in the 1932 temperature scale (Keesom 1932) it 
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is still not possible to fit all their results with a single value of 0, so that 
the anomalous peak still remains. It is believed that this behaviour must 
be attributed to experimental error in their results. 

Since nickel is ferromagnetic below 631°K, the heat capacity data include 
the contribution arising from the spontaneous magnetization. According 
to the Bloch spin wave theory, the heat capacity due to demagnetization 
for any cubic structure satisfies the equation (Mott and Jones 1936). 


C=1-18R 


Og—-o 
a (2) 
% 

where o is the magnetization at temperature 7' (7'<@) and oy is the satura- 

tion magnetization. The magnetization at low temperatures is given by 

the equation 
o=0,)(1—AT?) o> ee 
the constant A for nickel being found experimentally to be 

8-9 x 10-6 deg-?/2 (Fallot 1936). Combining (2) and (3), the heat capacity 

of nickel due to demagnetization is found to be 

C=8-8x 10-5 T%/2 joule mole-1 deg“ 

so that the contribution to the total heat capacity is about 2% at 2°K, 

rising to 3°4, at 4°K. The accuracy of the present experiment is insufficient 

to distinguish this effect from slight alterations in the values of y and @. 

If these parameters are corrected for magnetization effects it is found that 

y=7-05 x 10-3 joule mole—! deg-* 
0=468°K 

giving a @ versus 7’ relationship at low temperatures, which has the form 

shown by the dotted curve in fig. 4. While these corrected values are un- 

certain to a considerable extent because of doubt concerning the 
applicability of (2) to a real ferromagnet, it seems likely that the trend of 


the correction is as indicated, i.e. such as to make the @ versus 7’ curve 
have a more pronounced maximum at 7'=—0. 


§ 5. CONCLUSIONS 


Improved data for the heat capacity of chromium and nickel at liquid 
helium temperatures have been obtained by measurements on very high 
purity samples in a calorimeter employing a mechanical heat switch 
instead of exchange gas. The new results differ appreciably from those of 


previous observers and remove apparent anomalies in the @ versus 7’ 
curves for the two metals. 
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ABSTRACT 


It is shown that the well known dispersion formula of Drude is applicable 
to a dense medium even when the polarization fields associated with the 
different types of oscillators may all be different. The parameters 
involved in the formula are naturally the frequencies and the oscillator 
strengths of the oscillators in the medium, both of them as influenced 
by the mutual interactions between the oscillators, i.e. as influenced 
by the polarization fields. The observed data on dispersion will enable 
us to determine these parameters. In the Lorentz type of dispersion 
formula, on the other hand, the frequencies and the oscillator strengths 
involved are those of the isolated oscillators, which naturally are 
independent of the polarization field, and hence the effect of the polariza- 
tion field on dispersion has to be taken into account explicitly. Because 
of this, the Lorentz type of formula needs for its formulation explicit 
information regarding the polarization field, which cannot be supplied 
by the observed dispersion data alone. The reduction of the Drude 
formula to one of the Lorentz type, and the appropriate information 
regarding the polarization field that is needed for this purpose, are also 
discussed in the paper. 


§ 1. Tae DispeRsION FORMULAE OF DRUDE AND OF LORENTZ 


Tue well known dispersion formula of Lorentz for a dense medium is 


Kw—1l 4r2 N Fe? day ® a 


jets 3 2102-0) =: 3 &(Q2—a¥y’ say aes ed € 


where &’, is the dielectric constant of the medium for frequency w, N F 
is the number per unit volume, F; the oscillator strength, and p ; the 
reduced mass, of oscillators of frequency Q;. This formula is derived 
on the basis that the actual field acting on an oscillator in the medium is 
not merely the field in the medium but includes in addition the polariza- 
tion field, which is further taken to be 47/3 times the polarization per unit 


+ Communicated by the Authors. 
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volume. This introduces in the expression for K,—1 a multiplying 
factor [1+ (47/3)x.]=(K,+2)/3, where x, is the polarization of the 
medium per unit volume per unit field in the medium. 

On the other hand the Drude formula, namely 
, n N.f.e2 n . 
K,—l=4r5> ao Dies aa EAR oer £74) 


1 h(w?—o?*) t=10;_W 


which had been derived much earlier, did not naturally contemplate 
the presence of a polarization field. Hence the Drude formula is generally 
regarded as a special case of the Lorentz formula, applicable to media 
in which the polarization field is known to be zero. 

For convenience in discussion we shall refer to the quantities A, and 
a; briefly as the oscillator strengths of the oscillators concerned in the 
two formulae, though actually they refer to N, oscillators. 

Many years ago Herzfeld and Wolf (1925) showed, however, that 
formula (1) can always be reduced algebraically to (2), the w,’s and a,’s 
appearing in (2) being calculable functions of Q,’s and A,’s._ They derived 
(2) from (1) in the simple case when the number of frequencies involved 
in the formula is just two, but as they themselves pointed out, the 
reduction of (1) to (2) can be done also in the general case when the 
frequencies are more than two. 

Further since the number of parameters involved in the two formulae 
is the same, it should be possible conversely to reduce (2) to (1), though 
the actual reduction, as we shall see later in the present paper, is not 
quite so simple as the reduction of (1) to (2). 

From these considerations it follows that if the experimental data for 
the dispersion of any substance can be fitted with one of these formulae, 
they can be fitted equally well with the other formula too. The 
characteristic frequencies appearing in the two formulae will naturally 
be very different. From this circumstance it has been concluded by 
some of the later workers that by comparing the frequencies occurring 
in the two formulae with the observed frequencies of the medium it 
should be possible to decide between the two formulae : in other words 
to decide whether there is a polarization field of magnitude 47/3 times 
the polarization per unit volume, as contemplated in the derivation of 
the Lorentz formula, or it is zero as is taken to be implied by the Drude 


formula. 


§ 2. Tue DisperRston Data CAN GIVE NO INFORMATION REGARDING THE 
POLARIZATION FIELD 


In some earlier papers (Krishnan and Roy 1952, 1953) we have shown 
that such an appeal to observation, namely to find which set of frequencies 
deduced from the observed dispersion agrees with the observed resonance 
frequencies, whether it is the Q,’s or the w,’s, will not enable us to decide 
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on the existence or otherwise of a polarization field in the medium, and 
much less to determine its magnitude. Irrespective of whether a 
polarization field occurs or not, and irrespective of its magnitude when it 
occurs, the Drude formula (2) will always express correctly the observed 
dispersion, and the frequencies appearing in this formula will be the 
actual resonance frequencies of the medium. Since (2) can always 
be reduced algebraically to (1), again irrespective of the presence or not 
of an actual polarization field, the observed validity of (1) has by itself 
no physical significance, unless the magnitude of the polarization field 
can be verified independently, or the frequencies Q,, whose significance will 
be stated presently, are known from other considerations. 

This arises from the following circumstance. The effect of the 
polarization field, whatever may be its magnitude, on the dielectric 
constant, can be taken into account in two alternative ways, which are 
mathematically equivalent, and which lead to the Lorentz and the Drude 
formulae respectively : 


(1) The frequencies and the oscillator strengths may be taken to be 
those of the oscillators when isolated from one another, namely 2; and 
A, respectively, and the effect of the polarization field on the dielectric 
constant is then taken into account explicitly. 


(2) Alternatively, one may use the frequencies and the oscillator 
strengths, of the oscillators as influenced by the mutual interactions of 
the oscillators, i.e. as influenced by the polarization field, namely w,; and 
a;, in which case the whole effect of the polarization field on the dielectric 
constant is taken into account automatically. 

The difference between the two approaches lies in this: whereas w; 
and a; can be obtained directly from the dispersion data, 2; and A, 
cannot be. 


§ 3. ScopE OF THE PRESENT PAPER 


It was also emphasized in these papers that the polarization fields 
associated with the different types of oscillators present in the medium, 
as for example the electronic and the reststrahlen oscillators in the 
alkali halide crystals, may be widely different. It was further mentioned, 
without a statement of the proof, that the Drude formula (2) will hold 
even in the most general case when the polarization field acting on an 
oscillator of type i due to all the oscillators of type j in the medium is 
Pi; times x;, where x; is the contribution from oscillators of type j to the 
polarization per unit volume per unit field in the medium, and the Pig 8 
may all be different. Since this result is of some importance, and the proof 
is not obvious, the proof is given in the present paper. 

The conditions under which the Drude formula, which is generally 


applicable to any medium, can be reduced to one of the Lorentz type 
are also discussed in this paper. : 
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§ 4. THe Basic Expression FOR REFRACTIVITY AND ITS REDUCTION TO 
THE DruDE FoRMULA 
When the different p,,’s, which we may refer to as the polarization field. 


factors, may be different, the basic expression for the dielectric constant 
of the medium takes the form 


n 
K,—\=4ny=4r > xp (ule Uae a aon (CD) 
i=1 
where 
N ,F 7? Ww 
laren) eS 


Obviously the x’s, like K, would be functions of the incident frequency w, 
but for convenience in writing, we have dropped the subscripts w except 
in the case of K. 

We now proceed to demonstrate that the basic eqns. (3) and (4) can 
be reduced to formula (2) of Drude. 

Using the n equations similar to (4), and using also (3), one can eliminate 
all the x,’s and obtain therefrom 
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VN 
which is a polynomial of degree n in w?, which cannot be zero for all values 
of w?. 
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Now let P;; be the co-factor of the ijth element in the determinant 4. 
One then obtains from (5) 
n 
ioe cares >; P,/4 . . . . . . . . . (8) 
i,j=l 
In (8) the P,,’s are polynomials of degree n—1 in w*, whereas 4 is of 
degree n in w. Hence the right-hand side of (8) can be expanded into 
partial fractions, thus leading to the expression 


(= lesdary ss, lel eee (9 
K,—1=47x SP ley Te (9) 
where w,2, wo2, ... w,2 are the n roots of the equation A(w?)=0. w, and 


a, will be functions of the 2’s, the A’s and the p’s, and will be independent 


of w. 
Equation (9) can be immediately recognized as the Drude formula. 


$5. Tue Dispersion Data IN RELATION TO THE POLARIZATION FIELD 


This derivation emphasizes the basic nature of the Drude formula for 
dispersion. Far from being a special case of the Lorentz formula for 
zero polarization field, it is physically a more significant formula than 
Lorentz’s, even when there is a polarization field, which may be large and 
complicated. Just as for the isolated oscilators, their optical behaviour 
is determined uniquely by their characteristic frequencies 2; and their 
oscillator strengths A ;, so also the optical behaviour of a dense assemblage 
of such oscillators is determined by the corresponding quantities w; and 
a, characteristic of the medium. The parameters involving explicitly 
the polarization fields get eliminated. 

For the same reason the dispersion data by themselves cannot give 
any information regarding the polarization field. If the medium has 
n characteristic frequencies, the observed dispersion data, which obviously 
should fit into (9), can supply 2” parameters, namely the w,’s and the a,’s 
which would be just sufficient to determine the frequencies Q; and the 
oscillator strengths A; of the n oscillators, when all the p,;s are known. 

Now if the p,,;’s were all different their number would be n2, in which 
case, the p’s would be too numerous to be calculated from the 2n para- 
meters available from the dispersion data, even if all the Q,'s and the 
A,’s were known ; except in the trivial case when n=1. But actually 
the p,;'s are not quite so many, since in general p,,; should be equal to 
Pj; i which case the number of different p,,’s reduces to 4n(n+1). 
If the Q,s and the A,’s are known, then all the pi 8 Should be deter- 
minable from the dispersion data when n<3. With n—3 or even 
with n=2, it would be possible in practice to describe satisfactorily the 
dispersion over a fairly wide region of the spectrum. 

The result obtained just now may be of even deeper practical interest. 
Till now we have been concerning ourselves with the dispersion formula 
alone, i.e. the expression for K as a function of the frequency w, which 
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enables us directly to determine only the parameters w, and a;,, and in 
which the polarization fields do not appear explicitly. On the other hand 
the expression for K as a function of the density of the medium would 
involve the polarization field explicitly. In particular, for molecular 
substances like benzene, naphthalene and other organic compounds, 
in which the molecules retain their individuality in the condensed state, 
one can obtain the 2;’s and the A,’s from the dispersion in the vapour 
state, in which the polarization fields are negligible, and the corresponding 
w,;s and a,’s in the condensed phase from the dispersion data for this 
phase. A comparison of the two sets, referring respectively to the two 
phases, should enable us to determine all the p;,’s, provided we choose 
a spectral region in which 3 frequencies would be adequate to represent 
correctly the dispersion data. 

On the other hand the variation of dispersion with density, in the con- 
densed phase, which may be secured by varying the temperature, would 
not be so helpful, since the p’s also are functions of the density. 


§ 6. REDUCTION OF THE Basic FoRMULA TO ONE OF THE LORENTZ TYPE 


Having shown that the basic formulae of dispersion, namely (3) and 
(4), do reduce to Drude’s even when the p,,’s are all different, we may 
proceed to enquire whether in this general case, (3) and (4) can also be 
reduced to one of the Lorentz type, by which we mean a formula of the 


type 
o(K)=4r > aise ae es et (10) 


in which the Q,’s and the A,’s refer to the isolated oscillators, and have 
the same significance as in (1). The answer is that this can not be done 
in the general case when the p,,’s are all different. 

The problem is really one of finding an expression for K, i.e. for yx, 
which will not involve the component polarizabilities x; separately. 
One obvious method of doing this is to eliminate the y,’s between (3) and 
(4), and this procedure, as we have seen, yields the Drude formula, which 
involves the characteristic frequencies w,; of the medium, and not the 
frequencies 2; of the isolated oscillators, which we now wish to retain. 
On a closer examination of (4) and (3) one can see that the separate x;’s 
can be eliminated, without at the same time eliminating the ,’s also, 
only when the p,;’s have all the same value, say p. nit. 

When this is the case, the explicit expression for ¢(K) appearing in 
(10) becomes 


(11) 
or 
ee ee a, ae ee) 


where «=47/p—1, 
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$7. Repuction OF THE DRUDE FORMULA TO ONE OF THE 
LORENTZ TYPE 


We have shown that the basic formulae of dispersion, namely (3) and (4), 
can be reduced to the Drude formula (9) even when the p,,’s are all 
different, and that in the special case when they have all the same value p, 
(3) and (4) reduce to (12) also. Hence it follows when all the p;;’s have 
the same value, it should be possible directly to reduce (9) to (12), i.e. to 
obtain the constants Q; and A; in terms of w; and a; when «=47/p—1 
is known. The reduction is done in the following manner. 

Starting with (10), and remembering that «=4z/p—1, one obtains 
immediately 
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Now knowing w; and a; and knowing further that every p,; has the same 
value p=47/(2+1), we may put the genom on the right-hand 
side of (13)=0, and treating it as an aes in w’, obtain the n roots 
of the equation, namely M,?, M,?,...M,,?, and Pert obtain 


ee eh 


Since the M,’s are now known, the numerator on the right hand side can 


be put equal to p )C; JT (M,?—w?) where the C;,’s can be readily obtained. 
t g=1 
We thus obtain j#t 


in which the C';’s are obviously independent of w. 

Comparing (15) with (12) one can immediately see that if the basic 
assumption that every p,; has the same value, namely p, is correct, then 
every M, and C; appearing in (15) should be identical with every Q; 
and A ; respectively. 

But ‘the important point to emphasize is that the values of M, and C; 
thus obtained will represent correctly the frequencies and the oscillator 
strengths of the isolated oscillators only to the extent to which our 
knowledge of the polarization field is dependable. Since the dispersion 
data themselves do not supply any information regarding the polarization 
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field, one may choose any arbitrary «, and find the corresponding M , and 
C, which will make (15) algebraically identical with (9). Since the Drude 
formula (9) represents correctly the observed dispersion, (15) too, with 
the appropriate choice of the M,’s and C,’s corresponding to our choice 
of «, will represent equally well the observed dispersion, and this will 
be so irrespective of the actual values of the different p,,’s ; but then the 


M's and C,’s that appear in the formula will not be the actual Q,’s and 
the A,’s of the isolated oscillators. 
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XCV. Some Observations on the Solid State of Argon 
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SUMMARY 


A polycrystalline specimen of solid argon has been grown using a 
modification of Bridgman’s method; grains of up to 4mm diameter 
have been produced. A technique was used to etch the specimen 
thermally, and the results are described in this paper. 


In order to conduct experiments on the plastic properties of solid argon, 
some preparatory observations on the solidification of argon have been 
made. Polycrystalline specimens with a grain size of about 4 mim have 
been obtained by solidification from the melt, and thermally etched 
superficial markings observed. 

The purity of the gas used in these experiments was quoted by the 
suppliers as 99-93° , the main impurity being nitrogen. Initial experi- 
ments were made by reducing the temperature of a bulb containing about 
1 ml of liquid argon to just below its triple point (83-8°K). It was seen 
that argon formed ‘ vapour snakes’ of the kind observed in several 
compounds by Verschingel and Schiff (1954). The conditions given by 
them as necessary for vapour snake formation (viz. high triple point 
pressure, and a high value for the ratio between the latent heats of 
evaporation and fusion) are satisfied by argon, for which the triple point 
pressure is 516 mm of mercury, and the latent heats of evaporation and 
fusion are 1557 and 280 cals/g mole respectively (Guggenheim 1945). No 
special care was necessary for the snakes to form, the temperature of the 
walls above the argon surface not being higher than that of the lower walls. 
When only the lower part of a glass tube having two bulbs joined by a 
narrow vertical neck was cooled, the argon froze inwards from the walls 
until the neck was completely blocked by the solid. The contraction 
resulting from further cooling then caused vapour snakes to grow in the 
lower bulb ; a free solid-vapour interface is therefore not necessary for 
their formation. 

Transparent polycrystalline specimens were obtained by progressive 
freezing from the melt, as in Bridgman’s method. The argon tube was 
slowly lowered into liquid oxygen at about 70°K, and a stream of argon 
gas into the tube was maintained so that the solid grew at the rate of about 
a a EE ae ce 
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1! mm/minute. The solid was covered by a thin layer of liquid, the top 
of which was heated by condensation of the gas, so that the temperature 
gradient in the liquid was about 5deg/mm. This method of growth 
resulted in the production of polycrystalline specimens having mean grain 
diameters of about 4mm, the diameter also of the central part of the 
argon tube. Faster rates of growth caused smaller grain sizes, but 
subsequent spontaneous grain growth produced specimens similar to those 
resulting from slower rates. : 

It was found that the solid argon could readily be freed from the glass 
walls of the surrounding tube by immersing the tube in unstirred liquid 
oxygen at constant temperature, with the oxygen level above that of the 
argon. Sublimation of the solid next to the walls then took place, vapour 
condensing on the cooler parts of the walls near the oxygen surface, and 
the solid argon specimen becoming detached progressively, from the top 
downwards, from the containing walls. The freshly exposed solid-vapour 
surfaces rapidly showed thermally etched markings, which were visible 
when viewed by transmitted light, using a narrow aperture at the objective 
of the viewing microscope. 

The later specimens were grown in tubes having a narrow neck of 
rectangular (6 mm x 2 mm) cross section between bulbous ends, and the 
lower tail of the enclosing dewar vessel was of square section. By this 
arrangement the astigmatic effects of cylindrical walls were eliminated, 
and the argon surface could be studied more easily. The plate shows two 
photographs of typical solid argon surfaces. It should be noted that 
etched markings on both front and back surfaces are visible, particularly 
in fig. 1 (Plate). 

The observations of Chalmers et al. (1948) on the thermal etching of 
silver suggest that the strong markings to be seen in the photo- 
graphs are due to grain boundary grooves. That this is true in this case 
is supported by the observation of lines traversing the surface between the 
stronger markings, as in fig. 2 (Plate). These indicate that the strong lines 
enclose regions having a unique crystallographic axis, and we may there- 
therefore assume that these markings are in fact due to grain boundary 
grooves. Moreover, the fact that the grain boundaries do not show many 
kinks indicates that most of them are visible, as junctions with other grain 
boundaries usually cause sharp angles. This deduction is confirmed by 
the observation that the lighter surface markings in general terminate 
only at the visible grain boundaries. 

Grain boundary grooves are due to the excess free energy resulting from 
local atomic disorder, which gives rise to forces analogous to surface 
tensions (see Chalmers et al. 1948). The shape of any groove thus 
depends on the relative free energies of the grain surfaces in contact with 
each other and with the adjoining atmosphere. A few boundaries may - 
be seen to terminate at no other visible boundary, and fig. 2 shows 
surface striations bounded on one side by a line (AA) which is not itself 
visible. These observations indicate that the grains may in some cases 
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be twinned, the grooves formed at twin boundaries being shallower (and 
therefore not so clearly visible) because of the low free energy of the twin 
boundary. oe, 

Chalmers et al. observed the formation on silver of striations on the 
surfaces of the grains themselves, and showed that they were oriented along 
simple crystallographic planes. Such markings are due to the decrease in 
surface free energy sometimes made possible by the formation of steps, 
the faces of which are simple crystallographic planes. Similar striations 
are visible on argon surfaces, as shown particularly in fig. 2, but it may be 
noted that the lines are not quite straight. Figure 2 shows part of the 
surface of a specimen which had not been strained by application of an 
external load, but may have been subjected to considerable strains due to 
temperature changes in the specimen. These markings are not however 
slip lines, since they do not appear until the argon surface has been free 
for about twenty minutes. 

The conditions necessary for the visual observation of the surface 
markings should be noted, as they give some indication of the scale of the 
etching. Light from the illuminating source is deflected by the curved 
sides of the surface grooves, so that some rays which would have entered 
the microscope aperture are prevented from doing so. The grooves there- 
fore show as dark lines on a bright background, the contrast being im- 
proved by using a narrow slit source and a narrow aperture. The geo- 
metrical conditions necessary for surface steps to be visible are much less 
likely to be satisfied, the sides of the steps being plane. It is not surprising, 
therefore, that surface striations are visible on only a small number of 
grains. Further, for appreciable deflections to be caused, the grooves or 
steps must be at least one wavelength of light deep. 

It is not probable that the abrupt termination of some grain boundary 
marks is caused by non-satisfaction of the necessary geometrical conditions, 
as such an explanation involves pairs of conjoining grain boundaries being 
invisible, or else sharp kinks in some boundaries. If a grain is twinned, 
however, the grain boundary groove at the edge of one twin may be visible 
whilst its continuation at the edge of the other twin is too shallow to be 
seen. In a few cases, surface striations have indicated twinning in a 
region with edges marked by faint parallel straight lines extending across 
a grain. 

Preliminary observations of the properties of solid argon under applied 
stresses have also been made. The ultimate tensile strength at 75°K was 
found to be about 50 g/mm?, and, rather surprisingly, the final fracture 
diameter after necking was seen to increase with rise of temperature. The 
motion of grain boundaries during straining and subsequent annealing and 
recrystallization has also been observed (Domb 1955). 

This note is presented to indicate that specimens of solid argon having 
a grain size of several millimetres may readily be grown from the melt, 
and their surface etched to give information about their structure and 
about the behaviour of grain boundaries under internal or external stresses, 
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The thickness of the specimens grown in these experiments was only two 
millimetres, and Beck (1954) has reviewed evidence that the largest grain 
size obtained by spontaneous grain growth is limited by the smallest 
dimension of the specimen. However, later work has suggested that 
impurities are probably the most important factor determining the 
maximum grain sizes obtained in these experiments. There would 
appear to be no fundamental difficulties in the growth of single crystals of 
argon. 
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Wirrets (1953) and Wittels and Sherrill (1954) measured the density 
changes produced in quartz by fast neutron irradiation in a pile. The 
density was found to decrease with increasing dose, and the density 
change attained a saturation value of about 15°, for doses in excess of 
2 102° neutrons/em?. However, the rate of density change did not 
decrease steadily, but increased at first (beyond cumulative doses of 
7x 1019), reached a value of about two to three times the initial value 
at a cumulative dose of 10x 1019, and decreased only at higher doses—see 
also Kinchin and Pease (1955). To explain this increase in the rate of 
density change Primak (1955) suggested that the disorder is due to thermal 
spikes-regions whose temperature had been temporarily increased by 
the energy from the slowing down of a fast particle—and that these 
regions would retain their thermal energy longer in the partly disordered 
crystal because of the reduction in the thermal conductivity known to 
be associated with such disordering (Berman ef al. 1950, Berman 1951). 
Primak considered that this would enhance the density change produced 
by each new thermal spike. 

Primak’s explanation does not seem tenable when it is considered that 
the thermal conductivity of unirradiated quartz at the temperatures 
within the pile would be only two or three times the thermal conductivity 
after an irradiation of 7 x 101°, and that in any case a less rapid quenching 
of the temperature of a spike due to a decrease in the thermal conductivity 
would be expected to give rise to less, rather than more, disorder frozen 
into the spike. 

A different explanation of this increase in the rate of density change 
will be given here, based on the formation of disordered regions within 
the crystal, which are under compression until the number of such regions 
is sufficiently large that the stresses in the crystal are relieved by plastic 
flow. 

From the changes in the low temperature thermal conductivity on 
irradiation it can be concluded (Klemens 1951) that the imperfections 
responsible for the additional thermal resistance above liquid hydrogen 
temperatures are isolated imperfections of atomic dimensions, whereas 
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the additional resistance at lower temperatures is due to imperfect 
regions of larger size, though fewer in number. This picture is consistent 
with the theory of radiation damage (Seitz 1949). An atom, knocked 
out by a fast neutron, causes further displacements while being slowed 
down and stopped. In the initial part of the slowing-down process, 
the displacements are well separated spatially, but a substantial part 
of the energy is lost at the end of the range; where many displacements 
are concentrated in a small region. 

With increasing neutron dose the thermal conductivity curve approaches 
in character the curve for fused silica (Berman 1951), which suggests that 
within one of these larger disordered regions the material is of vitreous 
structure, as is reasonable from the glass-forming tendency of quartz. 

A similar picture of the vitrification of a region within the crystal 
by the slowing down of a particle was proposed by Primak (1954), who 
suggested that the substances which can be vitrified by fusion can also 
be vitrified by irradiation. Quartz is indeed turned into an amorphous 
statet by heavy irradiation (Wittels and Sherrill 1954, Primak 1955). 

The vitreous material formed by heavy irradiation is about 15°% less 
dense than crystalline quartz. If a small vitreous inclusion is formed 
inside the crystal, it will be unable to expand completely (i.e. 15°%), 
because of the constriction of the surrounding crystalline material. 
Thus the vitreous material will be in a state of compression, while the 
surrounding material will be displaced outwards, being sheared without 
dilatation. The resulting volume increase of the crystal can be calculated 
in terms of the elastic constants of both materials : for reasonable values 
of their Poisson ratios, and if their Young’s moduli do not differ greatly, 
the expansion is between one-half and one-third of what the expansion 
would have been if the vitreous material could have expanded fully. 

Thus for low doses the volume of the crystal will increase by 5-8% 
of the volume of the vitrified regions: at the same time shear stresses 
will be set up in the crystalline regions. When an appreciable fraction 
of the volume has been vitrified, the shear stresses at some points will 
be large enough to start plastic flow, and the stress system will be partly 
relieved. The mechanism of this stress relief is necessarily complicated 
and obscure, but the only way of relieving the stresses is by an increase 
in the volume of the crystal. In this way expansion will occur which 
was earlier prevented by the constraints of the surroundings, and the rate 
of density change due to irradiation will increase. 
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+ But this state differs from the state of fused silica, whose density is about 
2°, lower. Presumably the vitreous state is not unique, but different degrees 
of order are possible. Berman (1951) found small differences in the thermal 
conductivity of different specimens of fused silica. It would be interesting 
to compare the thermal conductivity of fused silica in the unirradiated and 
the heavily irradiated state. The latter, being denser, should be more ordered, 
and should have a longer phonon mean free path. 
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Primak et al. (1953) have studied the density changes of an irradiated 
crystal on annealing. Annealing at high temperatures increased the 
density to its original (crystalline) value: this is presumably due to 
recrystallization. However, annealing below 600°c caused a further 
decrease of the density. This is in agreement with the present picture : 
some plastic flow takes place when the temperature is raised, so that 
some of the stresses in the crystal are relieved and the latent volume 
increase realized, but the temperature is not high enough to produce 
recrystallization. 

It appears from an examination of the curve of density change versus 
dosage of Wittels and Sherrill (1954) that about one quarter of the 
material has been vitrified at a cumulative dose of 5x 101° neutrons/cm?, 
the density change being 2%. It is reasonable that the rate of density 
change should increase beyond this point, since the joining up of vitreous 
inclusions and the consequent formation of systems of high stress extending 
over large distances should be quite frequent. At twice the dose all 
single crystal x-ray reflections disappear. 

The fraction of material vitrified can also be estimated from Berman’s 
measurements of the thermal conductivity at low temperatures, and using 
the dosage values for these measurements deduced by Kinchin and Pease 
(1955), it is found that the two methods of estimating the fraction vitrified 
agree to within a factor 2, the thermal conduction values giving a lower 
fraction of vitritied material. The discrepancy could be due to a difference 
in the neutron flux calibration. 

According to the present model there is no direct relation between 
changes in hydrostatic density and changes in the lattice parameters 
determined by x-ray diffraction. Vitreous inclusions should produce 
very diffuse diffraction rings, as were indeed observed by Wittels. 
Diffraction spots and lines in the Laue and Debye-Sherrer patterns 
should become diffuse because of the inhomogeneous strain of the crystalline 
material, and long range order should progressively diminish. This is 
also in accordance with observations. However, apart from minor 
effects due to isolated point imperfections, there should be no net dilatation 
of the crystalline regions. The line broadening due to the inhomogeneous 
shear strain should be extremely skew, for in any volume element the 
compression along one principal axis should be twice the expansion in 
the other two directions. Thus the mean of a line should not be shifted, 
but the intensity maximum should be displaced in the sense of lattice 
expansion, to be compensated by a long tail of low intensity on the 
compression side. This could create the impression of overall lattice 
expansion and may explain at least a part of the lattice expansion 
reported by Wittels (1953) and by Johnson and Pease (1954). 
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ABSTRACT 

The existence of electron capture in the decay of 18F has been revealed 
by the detection of K-radiations of oxygen in a proportional tube spectro- 
meter containing a gaseous BF, source. The K/S+ ratio determined 
from the number of counts in the observed K-peak and the total number 
in the positron spectrum is 0-030, with an estimated accuracy of 6%. 
A scintillation spectrometer investigation showed the absence of any 
gamma rays with intensity greater than 0-5% of the positrons, suggesting 
that both K-capture and beta decay are ground state to ground state 
transitions. A comparison between theoretical and measured values of 
the K/B*+ ratio gives, for the ratio of the axial vector and tensor coupling 
constants, C,/Cp=(0-4+2)%. 


§ 1. IyrRopDUCTION 

THE ratio of the probabilities of K-capture and positron emission in 
an allowed transition is of interest for beta decay theory. Analyses of the 
measured shapes of allowed beta spectra have shown that Fierz inter- 
ference terms arising from the linear combination of different interactions 
are not large. If such terms are neglected, the probabilities A, for 
K-capture and A, for positron emission can be written 

G2 ; G2 
de= Za (Wot Wy GAME, = Sah(Z, Wo) LMP 
where G@ is the Fermi constant, Wy the maximum positron energy, W, 


the energy of the K electron, g, the large component of the wave function 
of the K electron at the nucleus, 


f(Z,Wo)= ie W(W2—1)12(W,— W)2F(Z, W)dW? 


with F(Z,W) the Fermi function, and WM is the nuclear matrix element. 
The K/f*+ ratio Ax/A,, does not involve the Fermi function constant and the 
matrix element, and is therefore particularly suitable for experimental 
verification. Moreover, agreement between experimental and theoretical 
values of the K/8* ratio can set a closer limit to the Fierz interference 
terms than that given by spectrum shape measurements. 

The ratios of K-capture to positron emission have been determined for 
a number of nuclei, and the results have been reviewed recently by 
Radvanyi (1955). However, the comparison with theory is not entirely 
satisfactory since the data concern nuclei with complex decay schemes, 
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not always well established. Further, indirect methods for measuring 
K-capture involving considerable corrections have usually been employed, 
and in most cases experimental errors are large. A K/8* determination 
for a simple ground state to ground state transition would provide a 
more reliable comparison with theory. In the case of a low Z nucleus 
such a measurement can be made directly using proportional counter 
gaseous source technique. With these considerations in mind, the 
present study of fluorine 18 was undertaken. 

Fluorine 18 is a positron emitter, of half-life 112 minutes (Blaser et al. 
1949), with a simple spectrum of allowed shape and maximum energy 
649 kev (Ruby and Richardson 1951). No gamma emission has been 
detected (Knox 1948). The transition is superallowed with log f7'=3-6. 
Fluorine 18 is an odd-odd N=Z nucleus with a spin of 1 (Moszkowski and 
Peaslee 1954), assigned on the basis of isotopic spin selection rules which 
indicate that the transition to '8O, with spin J=0 and isotopic spin 7'=1, 
involves a Gamow-—Teller interaction. 

Electron capture is expected, but it had not been observed before 
the present work, mainly because the absence of gamma emission prevents 
the application of the usual indirect methods for measuring K-capture. 
In 18F K-capture can only be detected by the K-radiation of oxygen 
subsequently emitted. This has an energy of 530 ev and, on account of 
the low fluorescence yield, 99:5°%, of the K-capture events lead to the 
emission of Auger electrons instead of K x-rays. <A proportional tube 
spectrometer can be used for the measurement of electrons in this low 
energy region, and this instrument is particularly suitable for the present 
experiment since the ratio of K-capture to positron emission can be found 
directly from the pulse height distribution. 


§2. SouRcE PREPARATION 


A radioactive gas to be measured in a proportional counter should not 
have a high attachment coefficient which might affect the performance of 
the counter, and of course should not react chemically with the wall 
material. For these reasons fluorine itself cannot be used, but preliminary 
experiments showed that boron trifluoride and methyl fluoride are more 
suitable. The former was employed in most of the measurements. 

The 18F was made by (y,n) reaction on fluorine using the gamma beam 
from the Glasgow 340 mev electron synchrotron. With a total output of 
2x 10° equivalent quanta per minute, the intensity of the effective part 
of the bremsstrahlung spectrum in the neighbourhood of 20 Mev is just 
sufficient to give a source of the required specific activity. 

Calcium fluoride, the target material, was irradiated for 4 hours and 
boron trifluoride was then prepared by heating the calcium fluoride with 
boric oxide and sulphuric acid. The gas was freed from condensable 
impurities by cooling with solid carbon dioxide. Since the target material 
was exposed close to the vacuum chamber of the synchrotron, the 
neutrons present produced an appreciable amount of *’A by (n, «) reaction 
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on calcium. This radioactive impurity was removed by freezing the 
boron trifluoride with liquid air and pumping out. 

The boron trifluoride prepared contained some silicon tetrafluoride 
which is undesirable in a counter because of its high electron attachment 
coefficient. However, the short half-life of 18F prevented the extensive 
purification required to remove the silicon tetrafluoride. The presence of 
this gas set a limit of about 4 mm of mercury to the partial pressure of the 
source that could be used in the counter. 


§3. MEASUREMENTS 


The proportional counter employed had a diameter of 10cm and a 
sensitive length of 58 cm, defined by field correcting tubes. A trace of 
gaseous source was introduced into the counter and methane to 10 cm of 
mercury, and argon to 3 atmos, was added. The total pressure was chosen 
so that the rise of the pulse spectrum at low energies, caused by escaping 
positrons that spend only part of their energy in the counter, would not 
be excessive. The pulse distribution was analysed with a multichannel 
kicksorter. 

Fig. 1 


COUNTS PER ENERGY INTERVAL 
n 
° 
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Kop X-RAY of OXYGEN 


«(2 ey ay 4 gah} ae ” ] j 
K-capture peak of 18F measured in a proportional counter with a gaseous 
source. 


i ape aate of the counter was adjusted using the 9-2 kev x-rays from a 

Ge sour ce for calibration, and the low energy part of the pulse spectrum 
was examined, A peak was detected in the region of 500 ev and is 
shown in fig. 1. To avoid the possibility of errors due to space charge the 
pOCley, of the observed peak was determined more accurately by using the 
2-8 kev K-radiation from an 37A source introduced into the counter. This 
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confirmed that the energy of the peak corresponds to the K absorption 
energy of oxygen. Such a peak can only be caused by K-capture since, 
in a counter with 47 geometry, any multiple event such as the emission of 
x-rays following ionization in positron decay would give an integrated 
pulse due to both components. The ratio of the counting rate in the peak 
to the counting rate in the continuous spectrum was found to remain 
constant as the source decayed, and the half-life agreed closely with that 
of '*F. Some long-lived activity due to impurities was apparent after 
several half-lives, but its intensity was low and it merely led to a slight 
increase in the continuous background. 

The relative intensity of K-capture was determined from the number 
of counts in the peak and the total number of counts in the spectrum. 
The peak is superimposed on a continuous positron pulse distribution, 
which rises at low energies and introduces some uncertainty in the position 
and shape of the line defining the base of the peak. Errors due to such 
uncertainties depend on the width of the peak, which, in low energy 
measurements, is relatively large because of fluctuations in the number of 
ion pairs produced. In the present case the width is increased even further 
by electron attachment in the source and by the effects of the high 
counting rate of large positron pulses. To reduce these errors a com- 
parison method was used when determining the area of the peak. After 
the K-capture measurement, when the source had decayed, ?7A was 
introduced into the counter and the shapes and relative areas of the 
230 ev L-peak and the 2:8 kev K-peak were measured. The base of the 
fluorine K-capture peak was established by comparison with these, and 
the number of K-capture pulses was thus determined. 

The high absorption of oxygen K x-rays in argon, and the small 
fluorescence yield, make any errors due to escape of photons from the 
counter completely negligible. The number of counts in the observed 
peak therefore gives the number of K-capture events that have occurred 
in the sensitive volume. 

The number of positrons was obtained from the total number of 
pulses recorded, taking into account background and the weak long- 
lived activity revealed by the decay curve. A correction, calculated to be 
4-6°%, was applied for positrons that enter the sensitive volume from the 
ends of the counter. 

Boron trifluoride is a reactive gas, so a separate experiment was carried 
out to measure any adsorption on the wall. The counter was filled with 
the source and, 1 hour later, was pumped out. The residual counting 
rate was then determined, and it was found that 23% of the ‘*F activity 
had been deposited on the brass wall of the counter. The decay curves 
taken with gaseous sources indicated that no loss of source occurred during 
the measurements, so it may be assumed that the adsorption takes place 
immediately after filling. Previous work has shown that such an adsorbed 
source would not contribute to the K-capture peak. The positron 
counting rate was corrected accordingly when calculating the K/pt 


ratio, 
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The results of a number of separate measurements give a mean value 
for the ratio of K-capture to positron emission of 
Ax /Ay=0-0380. 
The main uncertainty in the experiment arises in determining the number 
of pulses in the K-peak and, taking this into account, the estimated 
accuracy of the above result is 6%. 


§4. Dmcay SCHEME 

Although previous work (Knox 1948) has shown the absence of gamma 
rays above 100 kev with intensity greater than 10°, of the beta disintegra- 
tions, this does not rule out the possibility that the K-capture transition 
now observed might go to an excited state of 180. A more detailed search 
for gamma rays has therefore been made using a scintillation spectrometer. 
In view of the difficulty of detecting a weak photoelectric peak above the 
Compton distribution due to the annihilation radiation, the measurement 
was made by comparing the spectra obtained from !8F with that given by 
another pure positron source, in this case 14C. 


Fig. 2 
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Scintillation counter spectra of 18F and 110, 


A thin layer of irradiated calcium fluoride was placed in front of the 
Nal (Tl) crystal, which had a diameter of 4-2 cm and a height of 2-5 em, 
and was covered with a Perspex absorber. The pulse spectrum in the 
region between 30 kev and 700 kev was measured. The experiment was 
then repeated with a !!C source, prepared by (y, n) reaction from spectro- 
scopically pure graphite. Carbon 11 is a pure positron emitter with an 
end point energy of 990 kev and a half-life of 20 minutes. No gamma 
emission has been detected (Siegbahn and Peterson 1946). The spectra. 
obtained from both the !8F and 1C sources are given in fig. 2, normalized 
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to the same total number of counts. The close correspondence of the two 
spectra indicates that no gamma ray above 30 kev with intensity greater 
than 0-5 % of the positrons is emitted by either isotope. 

The energy region below 50 key was examined with the proportional 
counter and gaseous source. No peaks other than that due to the 
K-capture radiation were detected. It may therefore be assumed that the 


decay of '8F is simple, and the observed K-capture is a ground state to 
ground state transition. 


§ 5. DIscuSssION 


The theoretical value of the K/6* ratio for *F has been calculated using 
the table of Fermi functions by Rose (1955) with the screening correction 
factors of Reitz (1950). When the Fierz interference terms are neglected, 
the value Ry of the K/f+ ratio was found to be Ryp=0-0295, which is close 
to the measured value. 

By comparing the theoretical and experimental values of the K/B+ 
ratio it is possible to set a limit to the ratio C4/C,y of the axial vector and 
tensor coupling constants in the Gamow—Teller interaction. Studies of the 
shapes of allowed beta spectra, together with measurements of the angular 
correlation between electron and neutrino emission, indicate that C,4/Cy 
is not more than a few per cent. There are several estimates of the upper 
limit for this ratio: 0-02 (Konopinski and Langer 1953), 0-04 (Davidson 
and Peaslee 1953), and 0-20 (Winther and Kofoed-Hansen 1953). A de- 
termination of the K/S*+ ratio of ?*Na by Sherr and Miller (1954) gives 
C,4/Cyp=(—1+2)%, though there may be some doubt whether this decay 
is allowed. 

If the K/S+ ratio calculated for a small finite value of Cy/Cy has the 
value R, then it may be shown (see De Groot and Tolhoek 1950) that 

Cy R—Rk, 

Cy ~ 2Ro(1 + (W-?)) 
where < W~-") is the average of W~1 over the beta spectrum, W being the 
positron energy. For 1*F the value of (W~') is 0-69. Substituting the 
present values of the experimental and theoretical values of the K/p* 
ratio for R and R, gives C,4/Cyp=(0-442)%. 

This result is very similar to that obtained from the ?*Na measurements 
referred to above, and satisfactorily confirms conclusions drawn from the 
shapes of allowed beta spectra. 
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SUMMARY 


The angular properties of the reaction *Be («,n) !2C have been investigated 
in some detail in the range #,=0-4 to 1:3 Mev. Analysis of the results in 
terms of overlapping resonances has led to tentative assignments of spin 
and parity to certain levels in C0. The possibility that the reaction may 
proceed partly by a stripping mechanism has been considered ; it has 
been concluded that the results do not specifically demand such a 
mechanism. 


$1. INTRODUCTION 


THE suggestion by Madansky and Owen (1955) that the reaction 
®Be (a,n) *C may proceed by a stripping mechanism attaches renewed 
interest to some work done in this laboratory in 1951 and submitted in 
Ph.D. Theses by two of the authors (D.B.J. and G.A.J.) but not previously 
published. 

Work ante-dating our own included that of Halpern (1949) who 
measured the total neutron yield up to H,=5-3 Mev using «-particles 
from polonium, and Bradford and Bennett (1950) who, at #,—1-4 Mev, 
measured the 0° : 90° ratios of abundances of the two groups of neutrons 
to the ground state and first excited state of #C at 4-43 mev with #,,~6} 
and 2 Mev respectively. On the basis of these latter measurements, 
Haefner (1951) has assigned a spin of 5/2 to the relevant level in °C, but, 
as our own work shows, data at 0° and 90° are quite insufficient to make 
a reliable assignment. 

Post-dating our work, Talbott and Heydenburg (1953) have measured 
the y-ray yield up to H,=~3 Mev; they locate a well-pronounced 
resonance at H,=1-9 mev followed by a weaker, broader resonance at H,= 
2-6 Mev and a rising background. Bennett ef al. (1953, only in abstract) 
report the existence of a narrow resonance at H,—0-61 Mev and a broad 
one at H,—0-53 Mev. Tanner (1955) has measured angular distributions 
of the y-radiation at energies between H,—1:3 and 1-9 Mev and has 
estimated the total cross section at the 1-9 Mev resonance peak. 
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In this paper we present the excitation function of both groups of 
neutrons and of the y-radiation in the range H,=0-4 to 1-3 Mev, the angular 
distributions of the ground state neutrons at a number of energies in. this 
range, also, at fewer energies, the angular distributions of the low energy 
group, of the y-radiation and the n-y correlation. oie 

Attempts were made in the Theses to describe quantitatively in terms 
of three overlapping levels the angular distributions of the ground state 
neutrons. The qualitative success achieved did not encourage further 
application to the other distributions and correlations. Only the qualita- 
tive conclusions from these analyses are given here since in any case 
what was considered to be a single resonance at H,~0-6 Mev has been 
resolved into two (Bennett e¢ al. 1953) so that at least four overlapping 
states would be required to give a complete picture. When the analysis 
of angular distribution patterns demands such complexity as this there 
can be little profit in pursuing it further. 


§ 2. MEASUREMENTS AND RESULTS 


The neutron detector used in these measurements was a spherical 
ionization chamber of volume 500 cm‘, filled with deuterium to a pressure 
of 40 atmospheres. The chamber was constructed by G. H. Stafford 
and a similar one filled with hydrogen has been described by him (Stafford 
1948). A typical ‘spectrum’ from the reaction *Be («,n) C (thin target) 
is shown in fig. 1; two groups of neutrons are discernible with approxi- 
mate energies of cut off of 64 and 2 Mev respectively, and of relative 
intensities, after allowing for corrections listed later, of roughly 1: 3. 

For measuring the angular distributions of the neutron groups, the 
chamber was mounted on a rotatable framework centred upon the target 
where it subtended of solid angle of approximately 0-1 steradian. The 
yield was monitored by measuring the y-flux through a brass-walled 
Geiger counter at a fixed position—roughly forward. The target, 
an evaporated layer of beryllium on 0-010 in. copper sheet, was soldered 
to 1/32in. brass plate which formed the end of the target tube. The 
target thickness had previously been estimated from a comparison of the 
yield, from this target and from a thick target, of y-radiation from the 
reaction *Be (p, y) °B. Knowledge of the width of the broad resonance at 
L,=1 Mev (’=94 kev Fowler and Lauritsen 1949) gave a target thickness 
of 6 kev for protons of 1 Mev or approximately 40 kev for «-particles of 
1 Mev (Allison and Warshaw 1953). 

In measuring the angular distributions of the ground-state neutrons, 
pulses due to deuterons recoiling from neutrons of the low energy group 
were biased out. To detect this low energy group, a lower bias point was 
chosen in addition (just above noise due to build up of y-ray pulses etc.) 
and from the shape of the ‘ spectrum ’—see fig. 1—a suitable correction 
was made for the presence of the high energy group. Since this correction 
was only about 25%, of the total count errors in it are of little significance, 
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The angular distributions of the two groups presented in figs. 2 and 3 
have been corrected both for motion of the centre of gravity of the 
reacting system and for variation of detecting efficiency with neutron 
energy, which depends on the angle of observation. No correction for 
finite solid angle of detection has been made, because in general the 
statistical accuracy does not justify making this rather complicated 
correction. 


Fig. 1 
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The deuterium recoil spectrum of the neutrons from Be (a, n) ’C at #,=1-12 Mev. 


The detecting efficiency varies with neutron energy on four counts : 
(1) the n-d scattering cross section varies with the neutron energy ; 
(2) the deuteron angular distribution, and therefore the shape of the 
pulse spectrum, is a function of neutron energy ; (3) the discriminator bias 
cuts out a smaller fraction of the pulses at higher neutrons energies ; and | 
(4) the wall effect is larger for higher neutron energies, since the mean 
probability that a recoil deuteron be stopped within the chamber decreases 
as the length of the recoil track increases. 

The nd scattering cross section is given by Adair (1950). The effect 
of (2) is sufficiently small to be ignored. Corrections for (3) and (4) are 
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complicated and tedious to make and will not be discussed in detail ; 
they are not independent. 

The over-all correction for the ground-state neutrons resulted in an 
increase of about 20°% in the backward direction falling to zero in the 
forward, while for the low energy group the corresponding range of 
increase was 100% to zero. Though the corrections are large in this 
latter case, they are derived from well understood processes and should be 
reliable. 

The integrated cross section for emission of the ground-state neutrons, 
fig. 4, has been derived from measurements made at 0° and 60° and the 
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observed angular distributions. Those for the low energy neutrons and 
y-radiation, fig. 5, have not been corrected for the angular distributions, 
but, since the measurements were made at 90° and the respective angular 
distributions are not very strong, the correction is in any case quite small. 
It will be noticed that the three excitation functions are somewhat similar 
in shape indicating that the branching ratio between transitions to the 
excited state and ground state of C does not change very much with 
bombarding energy in our range. The yield in fig. 4 is given in terms of 
the absolute partial cross section for production of high energy neutrons ; 
it is based upon estimates of yields of 1-3 x 10~™ and 9x 107!" n per « at 
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E,=0-65 and 1-3 Mev respectively and upon the estimated target thick- 
ness of 6 kev to protons of 1 Mev, or 1-8 x 1018 Be nuclei per cm?, using a 
value of the stopping power of beryllium for protons given by Madsen 
and Venkateswarli (1948) and Allison and Warshaw (1953). The yields 
in fig. 5 are given in arbitrary units, but, as has previously been stated, 
the low energy neutron group is about three times as intense as the high. 

A thick walled brass Geiger counter, similar to that used for monitoring, 
was used as a detector in measuring the angular distribution of the 
y-radiation. The corrections applied here were for purely experimental 
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reasons such as for eccentricity of the geometrical arrangement and for 
absorption of radiation in the target backing. Since this part of the 
experiment was carried out in much better geometry this latter correction 
was found to be considerable when observing in a direction nearly parallel 
to the surface of the target, though very small in other directions. No 
correction was made for motion of the emitting nucleus. The results are 
shown in fig. 6. 

The angular correlation between the low energy neutrons and the 
y-radiation was measured only in the plane perpendicular to the «-particle 
beam. Liquid scintillation counters were used to detect both types of 
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radiation 7, the y-radiation chiefly by production of Compton electrons, 
and the neutrons by production of knock-on protons, within the scintilla- 
tion material (a solution of terphenyl in toluene). A coincidence circuit 
giving a resolving time of 1-2 x 10-8 sec was employed; the complete 
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+ Nal (T1) crystals for y-detection were not available in the laboratory at the time. 
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apparatus has been described by James and Treacy (1951). The need for 
fast coincidence techniques here was rather different than that usually 
encountered namely to discriminate against a high background rate from 
competing or spurious reactions ; the reaction now under consideration 
is virtually free from background but has a poor yield necessitating the use 
of large volumes of scintillating material for efficient detection. The 
rather poor light collection from these detectors made it necessary to run 
well into photomultiplier noise to realize this high sensitivity. 

To eliminate coincidences due to scattering of a neutron or y-ray from 
one detector into the other, measurements were taken only in the 90° 
180° sector. This eliminated completely the effect of neutron scattering 
by hydrogen nuclei since neutrons scattered through angles greater than 
90° (which requires multiple scattering) would have low energies and take 
longer than the resolving time to travel between detectors. The photons 
scattered through angles greater than 90°, and the neutrons scattered 
from carbon nuclei in the scintillators were absorbed by putting an inch 
of lead and a few inches of paraffin wax between the detectors. This 
was, of course, not possible at and near the 180° position; nothing 
was done about the neutron scattering from carbon, the probability of its 
occurrence in this position being quite small, but the y-ray scattering was 
effectively eliminated by surrounding completely one of the detectors by 
1 cm of lead, even in the direction of the target. This reduced the coinci- 
dence rate slightly (without affecting the correlation), but absorbed 
strongly the back-scattered photons, of energy about 250 kev. 

The angular correlations at two energies are presented in fig. 7. No 
correction has been made for finite size of detectors, which subtended 
+15° at the target. 

It might be remarked here that this correlation could be measured more 
accurately and reliably using Nal(Tl) crystals, now available, and slow 
coincidence techniques since the y-radiation could be efficiently detected 
with energy discrimination eliminating scattering events. 


§ 3. Discussion 

Analysis of the angular distributions of the ground state neutrons into 
powers of cos 6 is given in the table. The fact that odd powers of cos 0 
occur requires interference between states of unlike parity. In addition it 
has been found that the combination of even powers cannot be obtained 
without introducing interference between states of like parity. Thus at 
least three states of the intermediate system are required to interpret the 
results in terms of resonance reaction +. An attempt was therefore 


made to fit the angular distributions and the excitation function by the 
overlap of three resonances. 


gear i ae 
_ _} This conclusion is also apparent from the excitation function which 
indicates a resonance at or near 600 kev, while other work at higher energies 
ante and. Heydenburg 1953) indicate resonances in the neighbourhood of 
Mev of which the ‘ tails ’ of at least two may be important at low energies. 
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After exhaustive search of all reasonable sets of spin and parity assign- 
ments it was found that the combination (4-++) (4—) and (7/2—) for the 
relevant states in C was most promising and that the presence of a 
state of characteristics (5/2+-or—) was very unlikely since the ‘ pure ’ 
angular distributions from such states differed so much from the observed 


Analysis of angular distributions of ground state neutrons into powers 


. of cosd 
eee ee oe ee ee ee ee 

Ey ce Ge Os CG? Of 
(kev) 

510 1 —0-10 +0-04 = — = 

700 1 | +0-85+0-11 | —0-03-+0-06 | --1-28+0-14 — 

820 1 +1-64+0-10 | —0-02+0-07 | —2-44+0-13 — 

930 1 +1-22+0-08 | —1:39+0-20 | —1-:97-+0-11 1-45 0-18 
1040 1 0:81+0-07 | —2-16+0-18 | —1-53+0-11 | +2-58+0-18 
1090 1 0-63+0-12 | —2-47+0-29 | —1-33+0-18 | +3-11+0-31 
1140 1 | +0-30+0-05 | —2-98+0.14 | —0:79+0-08 | +3-73+0-15 
1240 1 | +0-07+0-14 | —2-92+0-34 | —0-69+0-20 | +3-84-+0-35 
1340 1 | +0-00+0-12 | —2-66-+0-33 | —0-43+0-20 | +3-63-+-0-37 
1140 1 | +0-21+0-12 | —2-81+0-33 | —0-73+0-20 | +3-61+0-36 

(repeat) 


patterns. However, as has been mentioned previously, Haefner’s tenta- 
tive assignment of spin 5/2 to the state at H,—1-9 Mev cannot be considered 
well founded. The complete expression for the differential cross section 
for this combination of states is, with cos 6 abbreviated to c : 


3 
W(8)= 55 A%(13—10c?+ 45c4) + 2. D4(9—3602+ 141¢4— 98c*) 


a s 4/5AD cos (—8) (3 —75c2-+ 2854 245¢8) +4. B24 40? 


3. fl AD 
2 4 al 


cos (B—8) (3—30c?+ 35c*) 
i 5 
—BC cos (S—y) +6] n4¢ cos (4—+y) (8c—5c?)—2 wv 


cos (8—-y) (83e—5c?), 

where A and « denote the amplitude and phase associated with neutron 
emission from the (3—) state formed by d-wave «-particles, D and 6 refer 
to the same state but formed by g-wave «-particles, B and f refer to the 
(4—) state, which can be formed only by d-wave «-particles, and C and y 
refer to the (4-++) state, which can be formed only by p-wave a-particles. 

The above combination of states can be made to yield the required 
experimental angular distributions at all energies if the relative amplitudes 
and phases are considered arbitrary. This is not the case, however, and 
on imposing the conditions that the yield curves for the individual states 
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and also the phases should have the usual restricted energy variation, 
dependent upon the location of the resonance, its width, and the barrier 
penetration function, great difficulty is experienced in fitting the data. 
The best fit is obtained when the (4+) state is assigned to the resonance 
near 600 kev, with its yield curve distorted by the barrier function to 
produce a considerable high energy tail. The contribution to the yield 
of the (3—) and (}—) states near and below this resonance is very small, 
but both rise steadily in intensity until at the high energy limit of our 
observations, they account for most of the yield. The analysis seemed 
to favour a sharper rise in intensity of the (}—) state giving a slight 
preference for assigning these characteristics to the resonance at 1-9 Mev. 

The full lines in fig. 2 are the theoretical angular distributions for this 
scheme and indicate that, though qualitatively they have the correct 
behaviour, much is still to be desired. The calculations were based upon 
a somewhat larger than normal nuclear radius for «-absorption (this 
point is discussed later) but for neutron emission the radius recommended 
by Feshbach and Weisskopf (1949) was used Tf. 

On adopting widths for the levels, the relative amplitudes and phases 
for all energies are defined by the resonance and Coulomb factors except 
for the relative importance of d- and g-wave for forming the ($—) state. 
Since relative barrier penetrability strongly favours d-waves, the g-wave 
contribution was neglected. 

The shape of the excitation function also leads to some difficulty since, 
using the conventional nuclear radius to calculate «-particle penetrabilities, 
it is not possible to produce the desired shape of the resonance near 
600 kev ; for such a slow rate of rise below resonance, the total width 
required is such that a minimum in the total yield would not occur. 
Neither does the yield rise sufficiently quickly at high energies. For these 
reasons a radius for °Be of approximately twice the conventional value 
was taken; subsequent estimation based upon the considerations of 
Scott (1954) on proton absorption indicate that such a choice may not be 
unreasonable. It is interesting to note in this connection that Thomas 
et al. (1949) give excitation functions for *Be (p, d) and (p, «) which seem to 
indicate no great increase in penetrability above a proton energy of 
about 350 kev. 

It is very likely that the resolution by Bennett et al. (1953) of what we 
have taken to be a single resonance near 600 kev into two, one broad, one 
sharp, may allow of a reduction in nuclear radius. The broad resonance 
will account for the slow rise at low energies, while the sharp resonance— 
blurred out by lack of resolution in our work—can be made to account for 
the slight dip in yield above the resonances. The rather slow rise at 
high energies may then arise from the higher contribution to the yield at 


| There is no necessity that the two radii should be the same. The analvsis 
seemed to depend uncritically on the radius for neutron emission so the con- 
ventional value was taken. 
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high energies of the broad resonance for which the resonance denominator 
works against the increase in penetrability. 

The intrusion of another state into the picture increases the possibility 
of fitting the angular distributions.+ No further attempt has been made 
by the authors to do this since the picture has now reached such a stage of 
complexity as to make it not worthwhile especially when the final con- 
clusions, namely the assignment of spins and parities to certain levels in 
*%C at excitation energies of about 11 Mev, are not likely to be of much 
value. It is clear however that our results are capable of being explained 
adequately in terms of the conventional resonance theory ; in particular 
no stripping mechanism of the type described by Madansky and Owen 
(1955) is demanded, but we cannot rule out the possibility that a small 
part of the intensity arises in this way. Although there is some tendency 
for the neutrons to be emitted backwards, this can easily be described as 
we have seen in terms of compound nucleus formation. Indeed the back- 
ward emission appears to depend too sharply on angle to be consistent 
with the stripping process but it is possible that stripping-compound 
nucleus interference may sharpen the backward rate of rise near 180°. 
The excitation function of Talbott and Heydenburg (1953) confirms and 
strengthens our conclusion since a non-resonant process such as stripping 
ean hardly produce what appears to be such a well formed resonance at 
1-9 Mev. 

Concerning the absolute yield at 1-9 Mev, from y-ray measurements by 
Tanner (1955), the assignment ({—), but not ($—), to the relevant state 
in 18C leads to an interpretation consistent with this measurement and 
also Dearnaley’s upper limit for o(«,«) quoted by Tanner. Taking 
200 kev for the width of this state, this latter measurement gives 30 kev 
as an upper limit for [’,, from which results a total neutron cross section 
at resonance of 0-18b. Allowing perhaps for 10-20% background 
intensity from other resonances gives a total neutron cross section at 
1:9 mev of 0-20-0-22b against Tanner’s measurement of 0:33 b to a 
factor of two. 

The assignment of (2) to this state and a partial width for «-particle 
emission of 30 kev at 1-9 Mev is not inconsistent with our measurement of 
the absolute cross section for emission of high energy neutrons at 1-3 Mev, 
but the attempt to correlate these data is very much dependent upon the 
penetrability function. Taking an effective nuclear radius of 6 x 10~ em 
the total neutron cross section from this resonance should decrease from 
0:18 b at 1:9 Mev to ~1-6 mb at H,=1-3 Mev. We find that the cross 
section for high energy neutrons at 1-3 Mev is 0:5 mb and for all neutrons 
approximately 2-0 mb. Since our analysis has indicated that this 
ee 2 ere he ad 7 et es a 

4 i r analysis was taken to be at resonance at Hy= 
a ie ene of the Se etion function. In fact the variation with 
energy of the angular distributions (see fig. 2) would suggest that this ponanes 
should be at a lower energy ; the location suggested by Bennett et a/. namely 
E,=530 kev would immediately lead to a better fit for the angular distributions. 
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resonance should account for roughly half the total yield at this energy, the 
agreement is, perhaps fortuitously, excellent. At H,=1-3 Mev, we 
deduce I°,~3 kev and the reduced width y,’~0-84 Mev cm or ~12% of 
the Wigner Limit. 

Assuming that the yield at #,~0-65 Mev arises chiefly from the (}+-) 
level, to which we have assigned a total width [=70 kev, we deduce that 
Ty~20 ev and y,2~0-5 Mev cm or ~7% of the Wigner Limit. In view 
of our complete lack of knowledge concerning the sharp resonance 
(except that, since we did not observe it, perhaps its intensity may be 
small), these latter figures are to be taken with reserve. Since both our 
estimates of y,? are rather larger than is normally encountered, it seems 
probable that we have not underestimated the effect of the Coulomb 
barrier despite the fact that a large effective radius has been taken. 

With no completely adequate account of the ground state transition, 
analysis of the other angular distributions seems profitless, but one point 
of interest arises from the n-y correlations. Rotation in the plane per- 
pendicular to the beam from 90° to 180° relative to the neutron axis 
corresponds to a rotation from 90° to ~60° relative to the impact axis 
under the experimental conditions. For a stripping reaction it is the 
impact axis which is of importance rather than the neutron axis. That 
the correlation should change so much in this range is additional evidence 
against the importance of the stripping mechanism. 


§ 4. CONCLUSIONS 


We consider that the experimental investigation of the reaction 
°Be (a, n)C up to #,=1-3 Mev reveals no clear indication of a stripping 
mechanism and that the main contribution to the yield in this region 
arises from four resonances: a broad resonance at H,—0-53 Mev corre- 
sponding to a level in the compound nucleus °C of characteristics ($+-) ; 
a sharp resonance at H,—0-61 Mev of unknown character; and two 
resonances at higher energies one of which is certainly that of H,=1-9 Mev, 
the other most probably that of = 2-6 Mev (see Talbott and Heydenburg 
1953) of characteristics (—) and (}—). In order to account for the 
cross section at H,—1-9 Mev we assign the characteristics (—) to the 
corresponding state in 13C, 
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XCIX. The Cold Work Introduced during the Fatigue of Copper 


By R. D. McCammon and H. M. ROSENBERG 
The Clarendon Laboratory, Oxford + 


‘(Received May 19, 1956] 


§ 1. INTRODUCTION 


Recent work (e.g. Wood 1955, Clarebrough et al. 1955) has suggested 
that the damage to the crystalline structure of a metal subjected to 
deformation depends upon the mode of cold working. In particular, the 
damage produced by fatiguing a metal is of a different nature from that 
produced by unidirectional strain. In the paper by Wood (1955) it is 
shown that the x-ray back reflection photograph of a polycrystalline 
copper specimen, fatigued in a compression—tension machine, differs 
from that obtained for a similar specimen strained in tension to the same 
degree of hardness. Furthermore, the temperature at which the damage 
anneals out, as determined by the x-ray photographs, is different in the 
two cases. The work of Clarebrough et al. (1955) showed that the energy 
which is stored in a metal when it is fatigued is released at a different 
temperature and at a different rate from that which is stored in a metal 
which has been strained in tension. The present note presents further 
evidence, from a different type of experiment in support of this view. 


§ 2. THe EXPERIMENTS 


The specimens were of polycrystalline high conductivity copper from 
I.C.1., annealed at 600°c for three hours in argon and furnace cooled. 
They were deformed in fatigue or tension by means of a Goodman’s 
vibration generator. By feeding a.c. or d.c. into this machine one obtains 
either an alternating (fatigue) or unidirectional stress. It has the 
advantage that the stress-strain curve of a specimen can be taken after it 
has been fatigued or strained without the necessity of handling the 
specimen between the two experiments. The degree of work hardening 
produced in a specimen after deformation is measured by taking a stress— 
strain curve at room temperature. 

Specimens were fatigued at various low temperatures at a given stress 
for 6-75 104 cycles (previous work showed that the work hardening — 
reached a saturation value after less than 104 cycles). It was found that 
the work hardening so produced decreased with decreasing temperature of 
fatiguing, even down to 4-:2°x—these particular experiments will be 
reported in detail in a future publication. The observations which are of 
interest in the present connection arose when a fatigue experiment was 
senna ment BO ee 
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done at 100°c (the specimen being actually immersed in boiling water) 
in order to see whether the work hardening introduced would follow the 
same trend as the earlier experiments and hence be greater than that 
introduced after fatiguing at room temperature. It was found that this 


was not the case, but that the work hardening was considerably less than 
the room temperature value. 


STRESS (DYNES /CM?x 10 7) 


4 8 12 
3 STRAIN % 


j 1 3 (a i t room temperature ; 
Stress-strain curves for copper specimens (a) fatigued a , | 
— (b) fatigued at room temperature and then heated to 100°C for 60 min. NS 

for 90 min x ; (c) fatigued at 100°c ; (d) virgin annealed specimen. 


It seemed improbable that this result could be ascribed to a reversal 
the fatigue characteristics of copper between room temperature an 
100°c. The most likely explanation appeared to be that some of the 
work hardening had annealed out at 100°c, although this was thought 
to be rather a low tempreature for any annealing to take place. To 
check this another two specimens were fatigued at the same loading a 
room temperature, heated to 100°c for 60 and 90 minutes respectively 
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and stress-strain curves then taken at room temperature. It was found 
that these curves almost coincided with the one obtained from the speci- 
men which had been fatigued at 100°c. These curves are shown in the 
figure, together with the stress-strain curve for a specimen which has been 
fatigued at room temperature at the same alternating stress but not 
annealed at 100°c. A stress-strain curve for the virgin annealed copper 
is shown for comparison. 

These experiments show that an appreciable amount of work hardening, 
although by no means all of it, anneals out at a comparatively low 
temperature which is far lower than the annealing temperature of the work 
hardening introduced by unidirectional strain. In order to confirm this, 
and also to check that the effect was not due to the apparatus, stress— 
strain tests were done on specimens which had previously been stretched 
so as to give approximately the same hardness as had been produced in 
the fatiguing experiments. It was found in this case that the stress— 
strain curves were the same whether the specimens were or were not heated 
to 100°c after the initial extension. 


§ 3. CONCLUSION 


It is thus apparent that the nature of the work hardening introduced 
by fatiguing is different from that produced by tensile straining. It is 
not as yet possible to say why the fatigue cold work partially anneals out 
at such a low temperature. At 100°c it is highly improbable that disloca- 
tions can themselves migrate. On the other hand it does seem that the 
dislocation configuration, which is produced by fatigue, is such that some 
of the ‘ pile-ups ’ involved can become unlocked at this temperature and 
hence produce the observed softening. 
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C. CORRESPONDENCE 


Energy Levels of A 
By I. J. van HEERDEN j and D. J. Prowse 
H. H. Wills Physics Laboratory, University of Bristol ¢ 


[Received May 29, 1956] 


VeERY few experiments have been performed to determine the energy 
levels of the *°A nucleus. They could be obtained from an analysis of the 
neutrons from the °K (d,n) 4°A reactions ; but the (d, p) reaction demands 
a supply of the rare isotope 9°9A. The 3’Cl («, p) 4°A reaction has been — 
investigated by Kranz and Watson (1953), but no levels in *°A were 
observed by them. It is also difficult to attack the problem through 
any form of resonance reaction. The easiest approach is to make experi- 
ments on inelastic scattering. Heitler et al. (1947) observed the scattering 
of 4-2 mev protons by argon, and found that the elastic scattering obeyed 
the Rutherford law to a high degree of approximation. They also found 
evidence for excited levels at 1:46 and 2:4mev. The first level has also 
been obtained from an examination of the y-radiation following K-capture 
of °K (e.g. Bell and Cassidy 1950). 

The elastic and inelastic scattering of 9-5 Mev protons by #°A has been 
studied by bombarding a gaseous argon target with protons from the 
beam of molecular hydrogen accelerated in the 60-in. Birmingham cyclo- 
tron. The angular distribution of the protons elastically scattered 
from *°A has already been given by Freemantle et al. (1954 b) ; in addition 
to the elastic group, an energy level was reported at 1-48-+-0-02 Mev. 

_ The experimental arrangement and procedure has already been 

described (e.g. Burrows ef al. 1951). Owing to the low cross section for 
inelastic scattering, it was necessary to work close to the scattering gap to 
obtain a sufficient density of tracks in the photographic plate. In this 
region of the plate, however, the locus of protons scattered through a given 
anlge is not parallel to the direction of the tracks themselves. 

The plate was scanned along the locus corresponding to the desired 
angle, and the projected track lengths were then corrected for * skewness ’. 
The true range of any particle in the emulsion was then deduced from its 
corrected track length and the angle of dip. To eliminate broadening of 
the groups due to single scattering, only those tracks which were nearly 
rectilinear and parallel to the correct direction were accepted. This 
procedure results in estimates of cross sections which are too low for the 
scattering process unless suitable corrections are applied, but it increases 


+ On leave from the National Physical Laboratory, Pretoria, 
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the resolution in energy. The tracks were measured on a Cooke, 
Troughton and Simms M 4000 microscope ; those > 100, with a ratchet 
device; those <100,, by means of an eyepiece scale. The loss of energy 
suffered by the scattered protons whilst traversing the gas was calculated 
from values of the stopping power of argon given by Brolley and Ribe 
(1955). The energy of the primary protons has been determined from the 
ranges of the elastically scattered protons, and the value thus obtained is 
9-51-0-01 mev. The range-energy relation used was that of Gibson 
et al. (1954). 

The number of background tracks (presumably due to particles 
scattered by the walls of the brass collimating tube), are at a minimum 
of angles of scattering of 115° (lab.). This is therefore a good angle to 
use in an experiment in which groups with a low cross section have to be 
resolved against a relatively high background. 


Fig. 1 
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Histogram of the tracks measured at 115° (lab.). 


A histogram of the tracks obtained at an angle of scattering 115° 
(L-system) is shown in fig. 1. It shows the presence of levels at 1-48 
0:02, 2:22+0-04, 3-12+0-03, 3-80+0-03, 4:50+0-05 and 4-98-+-0-05 mev. 
The errors quoted combine the calculated standard deviations in the 
measured ranges and beam energy, and allow for straggling effects. There 
is also a possible group at 2:66 Mev. Further measurements at 80° and 
130° (L-system) confirm the presence of the levels deduced from observa- 
tions at 1 15°. The 2-4 mev level observed by Heitler et al. may have arisen 
from incomplete resolution of the 2-22 and possible 2-66 Mev ‘levels shown 
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in fig. 1. The group at 1-48 Mev is in very good agreement with the 
previous work ; the other levels appear to be new. _ It is possible that some 
of these may be doublets ; thus, there are indications that the 4:50 mev 
level is accompanied by another at 4-40 Mev. The observed width of the 
group at 3-80 Mev suggests that this level also is a doublet. 

In view of the appearance of a large number of levels, the question of 
the purity of the argon used is of great importance. The argon gas was 
of normal commercial purity. Fortunately, the presence of impurities 
can be detected experimentally by observations on the groups scattered 


Fig, 2 
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Histogram of the tracks >300p in length, measured at 80° (lab.) showing the 
presence of impurity groups at ~380,. 


elastically by other nuclei. At all angles above about 40°, such groups 
from any oxygen or nitrogen atoms will be resolved from those from argon, 
but above 135°, these groups will begin to overlap with those due to in- 
elastic scattering by argon with the formation of the excited state at 
1-48 Mev. 

A histogram of the tracks obtained at'6=80" is shown in fig. 2. It is 
clear that there are impurity groups between those due to argon, but 
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it is not certain whether they are due to oxygen or nitrogen ; from the 
shape of the angular distributions it seems more probable that the principal 


impurity is nitrogen, as would be expected from on other grounds. The ~ 


absolute differential cross sections for the scattering of 10 Mev protons by 
oxygen, nitrogen and argon are known (Freemantle et al. 1953, 1954 a 
and 1954 b) and they enable us to deduce the percentage impurity in the 
argon employed. From measurements at various angles, it has been 
calculated that the maximum nitrogen impurity is 1-0--0-2%. From 
the known cross sections for the 2:3 and 3-9 mev levels of 14N, it then 
follows that less than 1 track due to the 2-3 level, and only 4 tracks due 
to the 3-9 level, would be expected to be present in the histogram shown 
in fig. 1. 

In addition, we have obtained indications of levels in *°A at energies 
>5 Mev (corresponding to ranges <100,). Further work is proceeding 
to confirm or deny the existence of these highly excited states. 
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Long Range Nuclear Interaction in Paramagnetic Resonance 


By J. E. Bennert and D. J. E. Incram 
University of Southampton t 


[Received April 12, 1956] 


Iv is well known that considerable information can be obtained from the 
hyperfine structure of paramagnetic resonance spectra, concerning both 
the nuclear parameters and the electronic state of the atom. In most 
cases so far investigated the hyperfine pattern has taken the form of 
(22-+1) lines, where J is the nuclear spin of the paramagnetic atom. 
Recent work (Owen and Stevens 1953) has shown, however, that ‘covalent’ 
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bonding between the central paramagnetic atom and its nearest neigh- 
bours may result in partial sharing of the magnetic electrons, and there is 
thus a contribution to the hyperfine pattern from these neighbouring 
nuclei. Such spectra have been analysed in detail for various chloro- 
iridate compounds (Griffiths and Owen 1954), and a similar interaction 
with surrounding halogen nuclei has also been found in some of the iron- 
group fluorides (Tinkham 1955). 

We have recently investigated a copper compound in which the halogen 
atoms are very far removed from the central paramagnetic atom, and yet 
have still obtained a marked contribution to the hyperfine pattern from 
them. The compound was copper tetra-phenyl porphin (Barnes and 
Dorough 1950), and the spectra of both the chlorinated and unchlorinated 
derivatives were studied. The phenyl groups are connected to the four 
carbon atoms round the edge of the central porphin ring system, and in 
the chlorinated derivative the outermost hydrogens of these four phenyl 
groups are replaced by chlorine atoms. The structural formula for the 
chlorinated derivative is shown in fig. 1, and it can be seen that there are 
seven atoms between each chlorine and the central copper atom. 

The hyperfine pattern observed from the unchlorinated derivative had 
the four equally spaced components of the normal copper spectrum, and 
the separations and g values were very close to those obtained for copper 
phthalocyanine, which has a similar central structure (Bennett and 
Ingram 1955). The pattern corresponding to the parallel direction is 
shown in fig. 2 (a), which is reproduced on the plate. 

When the chlorinated derivative was investigated, however, a very 
marked change in the pattern was observed and it had the standard 
christmas tree ’ appearance previously found in the chloroiridate spectra 
(Griffiths and Owen 1954). The pattern corresponding to the parallel 
case is shown in fig. 2 (b) and nine lines can be clearly resolved with faint 
signs of others in the wings. Since there was no appreciable change in the 
g value of 2-18, this type of structure can only be explained by assuming 
an interaction with other nuclei, and as the addition of the chlorines was 
the only change in the derivative, such an interaction would appear to be 
with the chlorine nuclei. Nitrogen nuclei are the only others present in 
the molecule possessing magnetic moments, other than protons, and any 
interaction with these would have been seen in both the unchlorinated 
derivative and the phthalocyanines. 

It is hoped that larger single crystals of the chlorinated derivative will 
enable a more detailed analysis of the spectra to be obtained, but a 
_synthesis of the pattern of fig. 2 (b) can be made if it is assumed that the 
magnetic electron of the copper is interacting with two chlorines at any 
given time, and that the chlorine hyperfine splitting is half that of the 
copper compounds, i.e. Agy=0-025 cm~* Aq=0-012 cm". Four over- 
lapping sets of seven lines are then obtained, each with an intensity ratio 
of 1: 2: 3: 4: 3: 2: 1, and these four sets are shown added together in 
fig. 2(c) with the envelope of the integrated absorption from them all. 
This has a fairly close resemblance to the photograph of fig. 2 (6) and would 
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Structural formula of chlorinated copper tetra-phenyl porphin. 
Fig. 2 (c) 
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Synthesis of pattern for chlorinated derivative, assuming chlorine hyperfine 
interaction half that of copper. Groups of 7 lines centred on the four 
copper components are shown below. 


seem to confirm that the interaction is with the chlorine nuclei. Although 
there are four equivalent chlorine atoms in the molecule the simultaneous 
interaction with only two might be expected as the resonating structure 
of the porphin ring has two sets of conjugated bond systems which are 
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equally probable. It is possible that the interaction im the parallel 
direction is with clorine nuclei of adjacent molecules. The large splitting 
observed shows that even in this case an unpaired electron must be asso- 
ciated with the chlorine for an appreciable time. 

It would therefore seem that the magnetic electron associated with the 
central copper atom can move out to the edge of very large molecules via 
the 7-orbitals of the ring system and may interact with nuclei present 
there. Although such long range interaction has been found in organic 
free radicals where the electrons move in non-localized orbitals, it does 
not appear to have been observed before with normal paramagnetic atoms. 
The apparent ease with which such interaction is possible may account 
for the strong exchange effects observed in many copper salts (Bagguley 
and Griffiths 1950, Bleaney and Bowers 1952). 
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The Bombardment of Osmium by Nitrogen Ions 


By G. W. GREENLEES and L. G. Kuo 


Physics Department, University of Birmingham 
[ Received March 28, 1956 | 


§ 1. [nrRopDUCITON 


Reactions caused by bombardment with heavy ions from a 60 in. 
cyclotron (maximum energy ~130 Mev) are likely to involve large spin 
changes because of the high values of incident angular momentum in- 
volved. Heavy ion bombardments could therefore be a means of 
exciting high spin isomeric states (Greenlees and Souch 1955). This 
could occur via a direct excitation mechanism and it was therefore 
decided to bombard a target with heavy ions and examine the activities 
of the target atoms theniselves after chemical separation. Osmium was 
chosen for these experiments since shell model considerations indicate the | 
likelihood of the presence of low lying levels of high spin. 


§ 2. PROCEDURE 


The bombardments were made using the internal nitrogen 14 beam of 
the 60 in. Nuffield cyclotron. The beam had a continuous energy distribu- 
tion of maximum energy about 130 Mev with an intensity fall-off by a 
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factor of two for each 10 mev rise of energy above 50 Mev. The target 
consisted of spongy osmium covered by a 5, foil of aluminium. This 
target was placed in a screened probe similar to that described by Chackett 
etal. (1954). A 5 aluminium screening foil was used. Beam currents of 
approximately 510-7. were available and bombardments of about 
1 A hour used. 

The osmium activities were separated chemically from the active 
target. The target was fused with sodium peroxide and the osmium 
distilled off as the tetroxide into alkaline solution. The osmium was then 
obtained as the sulphide by precipitation with H,S. This separation 
took about 14 hours so that any short lived osmium activities would not 
be observed in the separated fraction. The chemical extraction efficiency 
was about 70%. 


§ 3. EXPERIMENTAL RESULTS AND DISCUSSION 


The y-ray spectrum of the separated osmium was observed using a 
shielded NaI (T1) scintillation counter. A peak was found at 63 kev 
corresponding to the Ka x-ray energy for osmium. The decay curve for 
this peak is given in the figure. This curve analyses into three half lives 
of 7-2+1 hours, 39+-2 hours and 26+-2 days. 


1000 


COUNT /10 MIN. 
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Analysis of decay curve of 63 kev y-r ak smi 

_, Analy J ‘V y-ray peak from osmium. 

A, Half life 26+2 days; B, Half life 39+2 hours : C, Half life 7-2+-1 hour. 
Counts/10 min. Time-days. 
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These three half lives can be associated with some weak activities 
detected by Chu (1950). Chu bombarded rhenium with « particles and 
examined the iridium activities using a G.M. counter technique. On 
separating osmium from the iridium samples he found weak activities of 
half life 9-5 min, 6 hours, ~35 hours and ~50 days. The last three can 
be identified with the present results. The 9-5 min activity would not 
be observed in the present work because of the time required for chemical 
separation. In Chu’s work the only mass numbers to be expected are 
186, 187, 188, 189 and 190. All isotopes of osmium with these mass 
numbers are stable and hence the x-rays observed in the present experi- 
ments cannot be due to K-capture or the internal conversion of y-rays 
following 8 decay. This x-ray peak is therefore associated with the 
internal conversion of y-rays from isomeric transitions in osmium. ‘Such 
transitions might be expected to be highly converted which would 
explain the absence of the corresponding y-ray peaks with the scintillation 
counter. 

One concludes, therefore, that there are present at least three isomeric 
transitions in Osmium isotopes with mass numbers 186-190. From Chu’s 
work it is reasonable to associate the 7-2-+1 hours half life with Os!%, 
the 39+-2 hours half life with Os1*” and the 26---2 day half life with Os1®*. 
The absence of any y-ray half lives due to known radioactive osmium 
isotopes is evidence that the observed isomeric states are formed, in the 
present experiments, by a direct excitation mechanism. 
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Periodicity in Ice-Flowers 


By F. C. Frank and A. KELLER 


[Received May 16, 1956] 


THE accompanying photographs show ice-flowers formed on the inner 
surface of the panes of a cold green-house during the recent severe frost, 
and were taken on the morning of February 4th. 

Presumably the pattern developed on the evening of January 31st, 
on which day (following a period of mild weather) the temperature fell 
from 0°c at midday to —10°c in the night. The air temperature remained 
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below freezing during the following three days, but there were some periods 
of sunshine followed by similarly low night temperatures. 

The feature of interest is the regular periodic marking along the * stem ° 
of each ‘ fern-leaf’. The period was about the same in all “ fern-leaves” 
on any one pane, but varied with position in the green-house from about 
5mm to about 30mm. The periodic dark spots are regions of weaker 
light-scattering in the direction of view. This peculiarity in the frost- 
flowers was first noticed at night, on the evening of February 3rd, when 
the ice surfaces were more sharply etched, and at first sight gave the 
impression of an ice layer of considerable, sinusoidally varying, thickness. 
Holding the light in a suitable position, the banding while most prominent 
along the stem was seen to spread almost in semi-circles across the full 
width of each ‘ fern-leaf’. In the morning, when the sun was shining 
and some evaporation had occurred, it could only be seen on the stems. - 
The periodic bands were not visible between crossed polaroids, which 
enabled one to see, however, that the crystal orientation fanned out in 
the same manner as the surface striations, which were probably parallel 
to the hexagonal axes of the ice crystals. 

It is possible that this periodicity arose from a rhythmic crystallization 
analogous to Liesegang rings. Perhaps more probably it corresponded 
to a slow twist about the hexagonal axis (at the rate of 60° per apparent 
period) in the crystal growth. Such twisting, though unexplained, is a 
common accompaniment, usually with a much smaller period, of fanning 
crystallization, as in spherulites (see references). 
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CI. REVIEWS OF BOOKS 


Recent Advances in Optics. By. E. H. Linroor. (Oxford: Clarendon Press.) 

[Pp. ix+286.] Price 50s. 

THIs monograph gives an account of researches carried out in some selected 
branches of Optical Theory within the last twenty years. The first chapter 
deals with the Optical Image, more especially with the Geometrical Theory of 
error balancing, with diffraction images in the presence of aberrations and with 
the theory of partial coherence. Then follow chapters on the theory of the 
Foucault test, the theory of the Schmidt Camera and finally plate-diagram 
analysis with some of its applications, especially to plate mirror systems used 
in Astronomical Photography. 

The standard of presentation is unusually high. In particular, the handling 
or approximations so important in the theory of aberrations could hardly be 
improved upon. It is to be regretted, however, that from the wealth of material 
that undoubtedly comes under the title of this book, only very few topics are 
in fact covered ; much more could have been included in the space available if 
the topics were treated in less detail. Several small mistakes were noted by 
the reviewer: In equations (2.21), (2.22), (2.36) and (2.37) a multiplicative factor 
4iserronously included. It is not correct, as mentioned on p. 73, that the theory 
of partial coherence was first investigated by Van Cittert in 1934. Earlier 
investigations are due to Von Laue (1907) and Berek (1926). The author 
follows (regretfully, in the reviewer’s opinion) some recent practice of defining 
the degree of coherence without any reference to a time average; this will 
make it difficult for the non-initiated reader to appreciate the physical signifi- 
cance of this quantity. 

Much can be learned from this book not only about the subject matter it 
covers, but also about the elegance and precision with which theoretical optics 
can be presented. EW. 


Annual Review of Nuclear Science, Volume 5,1955. Edited by J. G. BECKERLEY, 
M. D. Kamen and L. I. Sceitrr. (Published by Annual Reviews Inc., 
Stanford, California; London agents, Messrs. H. K. Lewis, Gower Street, 
W.C.1.) [488 pp.] Price £3, or $7.50 postpaid. 


Tuts valuable series always covers a broad field of activities, and this time 
there is less than usual for the pure physicist. Four of the fourteen articles are 
on biological subjects connected with radioactivity, and four are chemical. 
These consist of one on methods of chemical separation, two on the chemical 
effects of irradiation, and one on the mass spectrometer in industrial chemistry. 

The physical articles include an especially clear and critical review of our 
understanding of nuclear gamma-ray transitions, by Goldhaber and Weneser. 
Peaslee writes on nuclear reactions, Ford and Hill discuss the distribution of 
charge in heavy nuclei, and Fretter reviews recent work on cloud- and bubble- 
chambers. 

In another interesting article, Borst offers critical advice to those who are 
wondering on what kind of reactor they should put their money (or their 
uranium). In discussing safety precautions, he remarks that there is, especially 
in universities, “ a tendency to avoid the routine and take short cuts ; and it 
is this virtue of the research scientist, making him outstanding in his field, 
which can constitute a significant and real danger.” 

This volume contains an index to the first five volumes, J.M,C.S, 
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The Fundamentals of Electroacoustics. By F. A. FiscuEr. (Translated by 
S Enruicu and F. Porpss. (Interscience Publishers London.) 


Tuts book, as the author informs us, resulted from the need to instruct young 
physicists and electronic engineers in the fundamentals of acoustics in order 
that they can work independently in this field as soon as possible. No special 
knowledge of acoustics beyond what is taught in this country at pass degree 
standard is required. The treatment is in the main theoretical though the text 
is written with an eye to practical applications and a bibliography of experi- 
mental treatises is appended. 

The author has certainly succeeded in his aim, though perhaps one should 
add that the reader whose training has been in electrical engineering will find 
himself more at home than one who has had courses in physics. The electrical 
aspects of the subject are in fact emphasized to some neglect of the hydrody- 
hamical aspects, so that the student might remain unaware, for example, of the 
effects of viscosity in causing damping. 

The translation is clear and the mathematics well presented. E.G.R. 


Basic Statistical Concepts. By Jok Kmnnepy Apams. (McGraw-Hill, 1955.) 
[Pp. xvi+304.] Price 41s. 6d. 


Tuts book is written by a professor of psychology with the aim of explaining 
the basic concepts of statistics and, at the same time, the mathematical language 
used in the subject: for example, elementary calculus. It is intended for 
students with no previous calculus or statistics. The treatment throughout 
is on an abstract level: applications, to psychology or to any other field, are 
rare. The concepts of power of a test and interval estimation are introduced 
early in connection with simple artificial examples dealing with finite popula- 
tions. They are later used in more mature situations. It is a serious criticism 
that the important concepts of efficiency of an estimate and the likelihood 
function are not mentioned. There is a long appendix on the calculus including 
derivations of the common sampling distributions. The book is not likely to 
please the mathematician; it will probably confuse the practical worker 
trained in British schools. DS Vele 


Tonized Gases. By A. von Ener. (Clarendon Press.) [Pp. xi+281.] 
Price 42s. 


Tats book fills a wanted place as a text-book for advanced students of the 
conduction of electricity in gases. It contains a great deal of information well 
arranged, especially on the individual processes that are the basis of discharge 
phenomena. Rather over two thirds of the book is concerned with these pro- 
cesses, which are studied in considerable detail. The experimental results are 
shown in a large number of clear diagrams, experimental methods are dealt 
with more sketchily. On the theoretical side a rough order-of-magnitude 
approach is wisely given wherever possible, and the more refined results are 
usually only quoted. Except in an appendix no attempt is made to deal with 
methods involving a change in the distribution function of electrons in a gas. 
Most of this book is concerned with conditions in which mobility is a valid con- 
cept, and more might perhaps have been said of the conditions of high fields 
and low pressures when this is no longer true. Regarded as a text book for 
the honours student, the balance, apart from this, is about right, but for those 
interested in research or in practical applications more on the actual discharge 
and less on the individual processes would have given a more useful book, 
especially as the latter have been so well dealt with by Massey and Burhop. 
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The style is rather abrupt, suggesting lecture notes, and at times the statement 

is a little too compressed, but the author is to be congratulated on his ability 

to get at the essentials of a physical argument and to put them across concisely. 
GiEeL. 


Progress in Nuclear Energy, Series 1: Physics & Mathematics, Vol. 1, 1956. 
Edited by R. A. Cuarpre, D. J. Hucuns, D. J. LirtLer, and S. Horowtrrz. 
With forewords by Sir Jonn Cockcrort and V. F. Wutsskorr. Published 
by Pergamon Press, London. [Pp. x+398.] Price 70s. 


Tuts is the first volume of one of six series of collected review articles based on 
the 1955 Geneva Conference on the Peaceful Uses of Atomic Energy, and 
including the results of much informal discussion made possible by the 
conference. 
A very considerable addition to the literature of neutron physics is made in 
the present volume, which contains the following reviews : 
Summary of Data on the Cross Sections and Neutron Yields of 02%, 
U?33, Pu®®, by J. A. Harvey and J..E. Sanders. 
Resonance Structure of U8, 17235, Py239, by P. A. Egelstaff and D. J. 
Hughes. 
Theoretical Analysis of Neutron Resonances in Fissile Materials by H. A. 
ethe. 
Techniques for Measuring Elastic and Non-Elastic Neutron Cross Sections, 
by L. Cranberg, R. B. Day, L. Rosen, R. F. Taschek and M. Walt. 
The Cross Section of Xe1°, by S. Bernstein and E. C. Smith. 
Resonance Capture Integrals, by R. L. Macklin and H. 8. Pomerance. 
Delayed Neutrons, by G. R. Keepin. 
Homogeneous Critical Assemblies, by D. Callihan. 
Physics of Fast Reactors, by J. Codd, L. R. Shepherd and J. H. Tait. 
Heterogeneous Methods for Calculating Reactors, by S. M. Feinberg. 
Highly Enriched Intermediate and Thermal Assemblies, by H. Hurwitz 
and R. Ehrlich. 
This publication is the first to contain a summary in this field of the work 
of scientists throughout the world, and will do much to spread the spirit of the 
conference as well as the fruits of the information presented. Gad: 


Physics and Chemistry of the Barth, Vol. 1. (London: Pergamon Press, 1956.) 
[Pp. viii+ 317.] Price 55s. 
Tuts is the first of an annual series, each of which will contain reviews of geo- 
physical and geochemical topics. Volume I includes :— 
. The Origin of the Solar System, by Sir Harold Spencer-Jones. 
. Temperatures within the Earth, by J. Verhoogen. 
. Radioactive Methods of determining Geological Age, by L. H. Ahrens. 
Seismology and the Broad Structure of the Earth’s Interior, by K. E. Bullen. 
The Hydrodynamics of the Earth’s Core, by Raymond Hide. 
Investigations under Hydrothermal Conditions, by Rustum Roy and O. iv 
Tuttle. 
7. Geochemistry of the Halogens, by Carl W. Correns. 
8. Geochemistry in the U.S.8.R., by 8. I. Tomkeieff. 
and Name and Subject indices. 
The series should provide valuable reference books for workers in these fields 
and will also be useful to other scientists who wish to follow recent developments 
in the study of the physics and chemistry of the Earth, Bace B 
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} } } . (Johns Hopkins 
Numerical Mathematical Analysis. By J. B. ScaRBOROUGH (Joh 
; University Press, London: Cumberlege.) [Pp. xix+554.] ae: 
Introduction to Numerical Analysis. By F. B. Hipesranp. (London: 
McGraw-Hill.) [Pp.x-+51l.] Price 64s. 


Turse two books are of a very different character. The first is the third 
edition of a work which originally appeared in 1930, well before the develop- 
ment of automatic digital computers had given the present impetus to the study 
of numerical analysis. It was, for a long time, one of the best of the relatively 
few books on the subject available. The author made substantial revisions for 
the second edition which was published in 1950. The present edition is largely 
a corrected reprint of the second edition with a little extra material on errors 


in the evaluation of determinants and a new chapter on the numerical solution — 


of simultaneous linear equations. 
The book by Hildebrand covers interpolation, quadrature, numerical solution 


of differential equations, least-squares polynomial approximations, the numerical — 


solution of equations, and related topics. Although the book is technically 
an introduction to the subject in that it starts at the beginning, it is really a text- 
book for the advanced student. It might, as the author suggests, be put into 
the hands of beginners to be read in conjunction with a course of lectures, but 
the treatment is really too detailed for those who are meeting the subject- 


matter for the first time. But for anyone who has progressed beyond this stage, 


and has had some practical experience either with a desk calculating machine or 
with a digital computer, the book is quite excellent. It is perhaps of most 
value as a source of theoretical results relating to the formulae and methods of 
numerical analysis, rather than as a source of advice for:the practical man 
faced with a specific problem. As an example showing the level of treatment 
sustained, I might mention the author’s discussion of the convergence with 
increasing of the Newton-Cotes n-strip quadrature formula. He shows that 
if the integrand has a singularity in the complex plain too near the axis of 
integration a divergent result is obtained, whereas the repeated application of 
the trapezoidal rule, or of Simpson’s rule, to the same number of strips may give 
a convergent result. 

The book is very well printed and edited, but it would be worthy of a more 
detailed index. _, M. Vows 


Progress in Metal Physics. Vol. 6. Edited by Bruck Cuatmers and R. Kina. 
(London and New York: Permagon Press.) [Pp. vili+354.] Price 70s. 


Tex volume for 1956 of this series, though it appears rather late, contains as 
usual much of interest. There are articles on the effect of pressure on resistance 
and on filamentary growth, a fashionable subject owing to the high strength of 
filaments. Mehl contributes an authoritative article on the austenite—pearlite 
raction and Sully an interesting article on creep which shows how little agree- 
ment there is between different investigators on the role of grain-boundary 
movement and cell formation. Knacke and Stranski have an excellent article 
on the mechanism of evaporation describing the known facts and with much 
information about the migration of atoms on a surface. Finally Hirsch provides 
a thorough and detailed review of our present knowledge of mosaic structure. 


N. F. M. 


[The Editors do not hod themse.ves responsible for the views 
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Thermally-etched grain boundaries on the front (A) and back (B) surfaces of 
a specimen of solid argon. (x 75.) 


Fig. 2 


The surface of a solid argon specimen showing thermally-etched detail within 
the grains, including the boundary AA. (x 90.) 
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